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(1) 


ABSTRACT 


The  thesis  deals  with  procedures  used  in  solving  two  types 
of  problems.  The  fixed  end  vibrating  string  is  nonlinear  in  that  the 
stress  strain  relation  is  perturbed  from  the  usual  elastic  form  to  give 
a  nonlinear  equation  of  motion.  The  second  problem  deals  with  an  elastic 
string  which  is  moving  through  two  fixed  eyelets  and  the  initial  conditions 
are  perturbed  from  a  zero  state.  Expressions  for  the  amplitude  of 
transverse  oscillations  are  derived  and  considerable  attention  is  given 
to  the  period  of  time  for  which  this  solution  exists  and  is  valid. 

The  techniques  used  involve  separating  the  problem  into  two 
separate  problems  both  having  characteristic  parameters  as  independent 
variables.  The  solution  in  terms  of  characteristic  coordinates  is 
obtained  by  finding  the  Riemann  invariants  and  expressing  the  derivatives 
of  the  amplitude  of  oscillation  in  terms  of  these.  The  mapping  of  the 
plane  of  the  characteristic  coordinates  into  the  space-time  plane  is 
approximated  by  calculating  first  terms  of  power  series  expressions  in 
the  perturbation  factor  e  .  That  is,  power  series  solutions  are  assumed, 
and  the  coefficients  of  the  system  are  also  expressed  in  power  series, 
and  then  the  equations  are  integrated. 

A  breakdown  of  the  transformation  from  the  characteristic  to 
the  physical  planes  is  observed  in  that  the  Jacobian  becomes  zero.  This 
indicates  the  boundary  of  the  interval  of  existance  of  the  solution,  as 
well  as  the  points  where  some  of  the  derivatives  of  the  amplitude  become, 
unbounded.  The  nature  of  the  mappings  at  the  point  of  breakdown  are 


considered . 
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Numerical  calculations  are  performed  to  determine  the  paths 
of  the  characteristic  curves  and  thus  determine  the  time  and  nature  of 


breakdown  in  the  case  of  the  fixed  string  problem. 
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CHAPTER  I 


INTRODUCTION 


The  purpose  of  this  thesis  is  to  invesitgate  the  solution  of  a 
number  of  almost  linear  hyperbolic  partial  differential  equations  in  one 
dependent  and  two  independent  variables  such  as  might  describe  the  motion 
of  vibrating  strings.  Many  of  the  techniques  employed,  in  particular, 
perturbation  procedures,  have  been  developed  for  problems  in  hydrodynamics. 
The  basic  idea  is  to  approximate  a  given  problem  with  a  solvable  linear 
problem  and  then  to  perturb  this  solution  in  terms  of  a  small  parameter. 
Solutions  are  expressed  in  terms  of  power  series  in  this  parameter.  We 
shall  attempt  to  correlate  and  explain  various  approaches  to  such  problems 
in  this  thesis. 

Breakdown  of  the  method  of  solution,  similar  to  that  associated 
with  the  formation  of  shock  waves  in  hydrodynamics,  occurs  in  these  problems. 
We  study  this  both  analytically  and  numerically. 

The  first  problems  to  be  considered  deal  with  a  vibrating  string 
with  fixed  ends.  The  string  is  nonlinear  in  the  sense  that  the  linear 
stress-strain  relation, 


T  -  T  (1  +  Y  ) 
o  x 


/ 


» 
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does  not  hold,  where  as  usual  T  denotes  tension,  is  equilibrium 

tension,  Y  is  lateral  displacement,  x  is  the  longitudinal  coordinate, 
and  subscripts  denote  partial  differentiation. 

Instead  we  replace  the  relation  (1.1)  by  one  of  the  relations 

T  =  T  (1  +  Y  +  \  e  Y2)  ,  (1.2) 

o  x  2  x 

or 

T  =  T  (1  +  Y  +  ^  e  Y3)  ,  ^1*3^ 

o  x  3  x 

where  e  is  a  small  parameter.  Then  we  obtain  the  equation  of  motion  by 
substituting  for  T  in 


P  Yt(;  -  Tx  .  (1.4) 

where  p  is  the  density  which  we  take  as  constant,  and  t  is  time. 

Attention  has  recently  been  drawn  to  these  equations  as  a 
result  of  numerical  experiments  carried  out  by  E.  Fermi,  J.  Pasta  and 
S.  Ulam  [7]  in  the  investigation  of  standing  longitudinal  oscillations 
of  nonlinear  beaded  strings.  We  briefly  describe  the  report  [7]. 

Consider  a  string  consisting  of  L-l  identical  point  masses, 
evenly  distributed  along  the  string  from  x  =  0  to  x  =  II  .  We  assume 
the  masses  are  connected  to  each  other  by  a  weightless  string  and  the  i 
particle  is  situated  at 
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x. 

l 


in 

L 


i  =  1,2, . . . ,L-1  . 


th 

Let  y_^  be  the  transverse  displacement  of  the  i  particle. 
The  equation  of  motion  corresponding  to  (1.3)  or  (1.4)  is 


~r  =  (yi+i  +  Vi  -  2yi)  +  0[(yi+i  -  yi)2  -  (yi  -  yi-i)2] 

d  T 


or 


(1.5) 


3t‘ 


1  =  (yi+i +  yi-i  -  2yi}  +  6[(yi+i  -  V  -  (yr  yi-i}  ]  ’ 


respectively.  Here  a  and  the  coefficients  of  the  quadratic  and  cubic 

terms  in  the  force  between  neighbouring  particles,  are  assumed  small,  and 
t  =  ct  where  c  depends  on  the  spacing  between  the  particles.  An  analysis 
for  various  times  t  is  attempted  in  terms  of  a  finite  Fourier  series: 


L-l 


kill 


r.  =  /  A,  (t)  sm  T 

i  i  i  k  L 

k=l 


(1.6) 


To  determine  the  coefficients  A^(t)  it  is  simpler  to  extend 
the  string  periodically  to  [0,2n]  .  We  can  then  use  the  relations 


2L-1  .. 

v  .  km  .j'Ui 

)  sm  —r—  sm  =  < 

•  1  1*  k 

1=1 


L 

0 


j  =  k  ^  0  , 

•  u  n  (1-7) 

J  T  0  , 

j  =  k  ^  hL  n  =  0,1,...  , 


so  that 
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i=l 


(1.8) 
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In  terms  of  the  y_^  variables,  the  total  energy  E  is 


ET  2^1  +  *  ‘  *  +  ^2L-1^ 

t  (2y^  +  ...  +  2y^_1  -  2yiy2  -  ...  -  2y2L_iy2L_2)  . 


(1.9) 


In  terms  of  the  A^.  we  have, 


l  2L;1  : 2 


2L-1 


ET  ■  1  l  K +  2L  l  <  -^'f 
k=l  k=l 


(1.10) 


If  we  associate  the  term 

1  !2  ,  ,.2  .2  nk  „ 

2Akt\m  2L  (1-U) 

th 

with  the  k  mode  of  vibration,  this  should  indicate  the  extent  to  which 

t  In  i 

the  k  mode  is  present.  The  conjecture  of  F.P.U.  was  that  the 

energy  should  be  equipartitioned  to  the  L-l  modes.  But  an  unusual  recurrence 

pattern  appeared  in  the  total  energy  of  various  modes.  After  a  certain 

length  of  time  t^,  most  of  the  energy  appeared  in  the  initial  mode  as 

depicted  in  Figure  I. 


K 


ENERGY  IN  LOW  MODES  OF  NONLINEAR  STRING 

Figure  I 

^F.P.U.  refers  to  E.  Fermi,  J.  Pasta,  and  S.  Ulam  with  specific  reference 
to  their  paper  [7]. 
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In  an  attempt  to  explain  these  phenomena,  Zabusky  and  Kruskal 
[15],  [17],  investigated  the  continuous  analogue  of  the  first  of  the 
F.P.U.  problems.  The  equations  they  considered  are  given  by 

0  <  x  <  l 

t  >  0 

0  £  x  £  l  (1.12) 

t  _>  0 

In  their  approach,  to  be  considered  in  Chapter  II,  a  breakdown 

4 

of  the  solution  was  predicted  at  time  t,  =  - -  ,  which  it  appears  is 

b  ealT 

too  small  to  permit  an  examination  of  the  recurrence  phenomenon.  The 

solutions,  which  are  in  effect  series  expansions  in  e  ,  are  for  t  <  t, 

b 

in  agreement  with  the  F.P.U.  results. 

A  second  estimate  of  the  time  of  breakdown  was  arrived  at  by 
P.  Lax  [9],  by  deriving  an  ordinary  differential  equation  involving  the 
derivative  of  one  of  the  Riemann  invariants,  and  showing  that  this  quantity 
becomes  unbounded  at  t  . 

A  procedure  related  to  the  one  used  in  [17],  and  similar  to 
methods  introduced  by  C.  Lin  [13]  and  P.  Fox  [9],  is  employed  in  this 
thesis.  It  involves  a  mapping  from  an  intermediate  plane,  a  plane  of 
"characteristic  coordinates',  to  the  physical  x,t  plane.  First  terms 


Y.  =  (1  +  eY  )Y 

tt  X  XX 


Y (x, 0)  =  a  sin  -  , 


i 


Yt(x,0)  =  0  , 


Y(0,t)  =  YU,t)  =  0 


of  a  series  representation  of  that  transformation  indicate  the  occurrence 


of  the  breakdown.  A  solution  to  the  second  F.P.U.  problem; 


Ytt  ■  (1  +  eYx)Y 

Y  (0 ,  t )  =  YU,t) 
Y(x,0)  =  f(x)  , 


xx 


0  <_  x  _<  il 
t  >  0 


=  0 

-  o 


(1.13) 


for  e  sufficiently  small,  is  found  and  the  corresponding  time  of  break- 
2 

down  indicated  . 


The  methods  used  in  the  fixed  end  vibrating  string  problem  are 
then  applied  to  another  class  of  equations,  those  describing  the  "moving 
threadline"  problem.  Consider  an  elastic  string  or  tape  moving  with 
constant  velocity  through  fixed  eyelets  or  rollers.  If  the  line  is 
displaced  by  a  small  amount  it  should  begin  to  vibrate. 


VIBRATING  STRING  PASSING  THROUGH  FIXED  EYELETS 


Figure  II 


2 

Because  of  the  breakdown  of  the  transformation,  the  analysis  of  the  above 
equations  cannot  yeild  a  result  corresponding  to  the  recurrence  observed  in 
the  F.P.U.  report.  Kruskal  and  Zabuski  have  since  published  a  paper  dealing 
with  a  more  comolex  equation  whose  solution  does  correspond  to  the  results  of 
that  report.  It  appears  in  'Interaction  of  Solutions  in  a  Collisionless 
Plasma  and  the  Recurrence  of  Initial  States',  Zabuski  and  Kruskal,  Physical 
Review  Letters,  Vol.  15,  9  August,  1965. 
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It  is  of  interest  to  determine  both  the  motion  of  the  line  and  the  time 
when  breakdown  of  the  solution  might  be  expected.  Physically,  breakdown 
might  correspond  to  a  snapping  of  the  line.  We  consider  two  models 
describing  this  situation  as  developed  by  J.  Zaiser  [18]  and  presented  by 
Ames  [1].  These  are  given  by 


^7  Y  +  ctY  .  + 

4  xx  xt  tt 


XX 


1+Y' 


X 


I  (1.14) 


a 


4 (1+Y  ) 
x 


Y  + 
xx 


a 


(1+Y  ) 
x 


-  Y  +  Y 

2^  xt  tt 


XX 


1+Y' 


x 


II  (1.15) 


Y  (0,  t)  «  Y(l,t)  =  0  ,  t>0  , 


Qi 

where  -j  denotes  the  ratio  of  the  speed  of  the  string  to  the  speed  of 


sound  in  the  medium,  t  is  time,  Y  is  lateral  displacement,  and  x  is 
longitudinal  coordinate.  Equation  (1.15)  differs  from  (1.14)  in  that  the 
longitudinal  velocity  of  individual  elements  is  varied,  that  is  is 
replaced  by 


a 


2  (1+Y2)15 
x 


The  perturbation  parameter,  in  these  problems,  unlike  the  fixed 
end  problems,  does  not  occur  in  the  equation  of  motion;  we  shall  introduce 
it  into  the  initial  conditions  as  follows: 


Y(x, 0)  -  e  sin  IIx  ,  0  <  x  <  1  , 

Yt(x,0)  -  0  . 


(1.16) 


Extensive  numerical  results  are  available  for  this  problem  in  J.  Zaiser' s 


• , 
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thesis  [18].  For  values  of  t  for  which  the  calculations  were  carried 

out,  no  singularity  in  the  solution  or  its  derivatives  was  observed  in 
[18]. 


In  this  thesis,  expressions  for  the  first  terms  of  a  solution 
in  power  series  in  e  are  developed  and  it  is  shown  that  our  solution 
does  exhibit  breakdown. 


The  breakdown  is  in  fact  predicted  for  a  value  of 
by  Zaiser’s  calculations. 


t,  covered 
b 


We  deal  with  two  forms  of  pertubation:  first,  perturbation  of 
the  equations  of  the  system  and  second,  of  the  initial  conditions.  We 
show  in  Appendix  I  that  these  in  general  are  equivalent. 


We  now  include  a  brief  discussion  of  the  various  approaches 
available  to  us.  We  are  concerned  with  equations  of  the  form: 


aY  +  bY  +  cYt  =  0 

XX  xt  tt 


(1.17) 


in  the  half  plane  on  t  >  0  ,  where  Cauchy  data  is  given  on  t  =  0  , 


Y  (x  ,0)  =  f  (x) 

—°0  <  x  <  •+«> 

Yt(x,0)  =  g (x) 


(1.18) 


with  f(x)  continuously  differentiable.  The  coefficients  a,  b,  c  are 
functions  of  Y  ,  Y  only  and  perhaps  of  a  parameter  c  .  Such  a  condition 

X  t 

in  the  coefficients  opens  a  number  of  approaches  and  is  necessary  to  the 
ready  solvability  of  our  equations  by  the  techniques  we  shall  use. 


1  fane  bsqolsvsb  sib 
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Since  the  data  for  the  actual  problems  we  consider  is  given  over 
a  finite  space  interval  0  _<^  x  _<  1  ,  it  is  necessary  to  extend  our  data 
by  a  continuation  to  the  whole  x  axis. 

We  write  the  system  as  two  first  order  equations  by  introducing 
two  variables  u  =  Y  ,  v  =  Y  ,  and  writing  the  equivalent  matrix 

X  ti 

form  [4] } 

0 

C  i  0  ,  (1.19) 

where  f(x)  =  -3—  f(x) 

dx 

The  procedure  used  in  this  thesis  consists  in  writing  two  pairs 
of  differential  equations  for  u,v  and  x,t  with  parameters  £,p 
(characteristic  variables)  as  independent  variables.  One  might  consider 
the  use  of  a  number  of  other  techniques  in  the  solution  of  the  problems. 

Because  the  coefficients  are  functions  of  u  and  v  only, 

the  hodograph  transformation  would  seem  useful  in  order  to  obtain  a  linear 

system.  This  involves  writing  the  system  with  x,t  as  dependent  variables 

of  u,v.  But  the  mapping  of  (u,v)  into  the  (x,t)  plane  is  many  valued 

except  in  restricted  ranges,  which  further  complicates  the  assigning  of 
boundary  conditions.  Although  the  coefficients  are  functions  of  the 


* 
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independent  variables,  finding  a  solution  is  still  not  trivial. 

Zabusky  has  solved  the  linear  equation,  obtained  by  a  hodograph 
transformation,  by  the  use  of  Riemann* s  method  for  the  first  fixed  end 
problem  [15].  He  was  able  to  write  down  the  solution  which  is  locally 
valid  because  the  required  Riemann  function  is  known,  and  by  an  "unfolding 
process",  the  solution  is  extended  to  the  half  plane.  However  application 
of  this  procedure  is  limited  by  the  scarcity  of  known  Riemann  functions  [2], 


An  approach  similar  to  the  one  we  use  consists  in  attempting  a 
solution  in  the  form: 

1 


IU1 

= 

(U1 

1  +  £(u) 

\V  1 

\  V  J 

I  \  V  / 

+  .  . 


(1.20) 


If  we  can  write 


A(u,v)  =  Aq  +  sA^  +  e  A^  +  . . .  , 


(1.21) 


with  Aq  constant  we  obtain  a  sequence  of  problems 


(o) 


IT 

.  /  u  \ 

+  A  = 

o\v| 

’  t 

X 

C)<n>  +  *o G)(n>-- VjKI* 

t  x  3=1  J  x 


(n-j) 


n  >  1  , 


(1.22) 


with  initial  conditions  u^n^(x,o)  =  f^  ^ (x)  ,  v^  ^(x,o)  -  g^  ^ (x)  , 


CO 


where  f(x)  =  .  £  e^f^^(x)  and  g(x)  =  ,  E  e^g^^  (x)  .  Since  the 

1C“  O  K.  O 


left  hand  side  of  each  equation  is  linear  with  constant  coefficients, 
solutions  of  the  sequence  of  problems  may  be  carried  out  successively.  The 


/ 


' 


11  - 


difficulty  with  such  an  approach  in  problems  where  breakdown  occurs  is  that 
it  is  difficult  to  obtain  information  in  the  neighbourhood  of  singular 
points  because  the  solution  series  often  do  not  converge  there.  The 
advantage  of  the  method  we  use  is  that  the  series  obtained  converge  up 
to  and  beyond  the  breakdown  point  [9]. 


'  ' 


CHAPTER  II 


PERTURBATION  METHODS  OF  SOLUTION;  CHARACTERISTIC  VARIABLES 


Consider  the  system 


lit  +  A>Jx  =  0 


(2.1) 


where  A  =  A(^a)  is  a  2  by  2  matrix,  and  u  = 


is  a  two-dimensional 


vector.  The  suffixes  denote  partial  differentiation.  We  suppose  that  A 
has,  in  the  region  under  consideration,  two  real  and  distinct  eigenvalues 
P,  Q,  so  that  the  system  is  hyperbolic.  Let  nu  =  (nu,n_^)  ,  i  =  1,2 
be  the  real  left  eigenvectors  corresponding  to  P,  Q  respectively,  so 
that  we  have: 


m^A  =  Pm^ 


%  v 

m^A  =  Qm2 


(2.2) 


It  immediately  follows  from  equation  (2.1)  that 


'Xi  v 

m.u  +  Pm-U  =  0 
l^t  lvx 

V  >x, 

m^u  +  Om~u  =  0  , 
2  vt  •  2^x 


(2.3) 


or  equivalently  that 


m,  (u  +  Pu  )  +  n.,  (v.  +  Pv  )  =  0 
1  t  X  I  t  X 


m2(ut  +  Qux)  +  n2(vt  +  Qvx)  '  0  • 


(2.4) 
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Because  A  is  a  function  of  u,  v,  e  only,  so  are  P,  Q, 

m.  ,  n,  ,  i  =  1,2 
11 

We  now  introduce  the  characteristic  variables  £  =  £(x,y)  , 
n  =  p(x,y)  .  We  assume  a  suitably  smooth  solution  (u,v)  and  stipulate 
that  along  the  curve  £  =  constant  that 

x  -  Pt  =0  ,  (2.5a) 

n  n 

and  along  p  =  constant  that 

x  -  Qt  =  0  .  (2.5b) 

We  assume  that  we  can  do  so  in  such  a  manner  that  the  Jacobian  is  non¬ 
zero  : 


J 


3(x,t)  /  „ 

3(5, n)  9 


In  this  case,  the  equations  (2.4)  can  be  written  in  the  form: 


1  u 

+ 

n.  v 

l  n 

l  n 

+ 

n  v\ 

2  £ 

2  £ 

0 

0 


along  E,  =  constant 
along  p  =  constant  . 


(2.6) 


The  equations  (2.5),  with  the  properly  transformed  initial  conditions, 
describe  the  mapping  from  the  (x,t)  to  the  (£,q)  plane  while  (2.6) 
describes  that  between  the  (u,v)  and  (£,n)  planes. 

We  designate  the  characteristic  variables  by  setting  the 
parameters  £,q  equal  to  the  x  coordiante  of  the  point  at  which  the 
characteristic  curves  n  =  constant  and  £  =  constant  interesect  the 


♦ 


14  - 


initial  line,  t  =  0  .  In  our  work,  the  slopes  P,  0  of  the  characteristic 
lines  in  the  (x,t)  plane  are  non-zero.  Thus  the  point  of  intersection 
of  a  particular  line  with  the  x-axis  is  unique. 


CHARACTERISTIC  LINES  £  =  const. ,  p  =  const. 

IN  THE  (x , t )  PLANE 

Figure  III 

In  this  way,  the  initial  conditions, 

U(X,0)  =  f(x)  -co  <  x  <  °°  , 

(2.7) 

v(x,0)  =  g(x)  , 
may  be  rewritten  in  the  form: 

u (£ , £  )  =  f  (£) 

(2.8) 

v  (£ ,  £ )  =  g(£) 


00  >  X 
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x(£,£)  =  l 


(2.9) 


tU,0  =  0 


In  the  problems  dealt  with  here,  as  in  the  pacers  by  Fox  [9] 
and  Lin  [13],  the  coefficients  of  the  system  are  functions  of  u,  v 
and  e  only  so  that  (2.6)  can  be  solved  independently  of  (2.5).  Methods 
which  utilize  the  factor  e  to  perturb  from  a  simpler  state  generally 
are  powerful  and  useful,  but  we  use  the  Riemann  invariants  associated 
with  the  characteristic  lines  because  of  the  simplicity  and  directness 
of  the  method. 

To  obtain  the  Riemann  invariants,  let  us  suppose  that  we  can 
find  two  non-zero  functions 


hl  =  h]_(u’v)  »  h2  =  h2('U,V‘)  * 


(2.10) 


and  suppose  further  that  there  exist  functions 


r  =  r(u,v)  ,  s  =  s(u,v)  , 


(2.11) 


such  that 


r  =  h-m-  , 
u  11 


(2.12) 


If  this  is  possible,  the  system  (2.6)  becomes 
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r  =  r  u  +  r  v 

n  u  n  v  n 


=  0  along  £  =  const. 


S_  =  S  U,.  +  S  V. 

£  u  £  v  * 


=  0 


along  n  =  const. 


(2.13) 


which  shows  r  and  s  are  constant  along  lines  £  =  const.  and 
n  =  const,  respectively.  Thus  one  concludes  that 


r  =  r (u,v)  =  r  (£) 
s  =  s(u,v)  =  s(n) 


(2.14) 


and  we  can  calculate  r  and  s  for  points  in  the  plane  by  doing  so 
at  the  initial  line. 


CHARACTERISTIC  LINES  IN  THE 
(£,n)  PLANE 


Figure  IV 


, 
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Along  A  P  , 


r  =  r(^1)  =  r(u(^1,C1),  v^,^)) 

=  r(f  (C1) ,  gC^)) 

s  =  str^)  =  s(f(n1),  gt^))  • 


(2.15) 


Thus  r  and  s  are  known  for  P(£^,ri^)  and  in  general  we  have: 

u(5,n)  =  u(r(£) ,  s(n)) 
v(S,n)  =  v(r (C) ,  s(n))  • 


(2.16) 


The  feasibility  of  this  method,  and  its  applicability  to  our 
problems  depends  on  the  fact  that  the  coefficients  are  functions  of 
u,  v,  e  only  and  that  the  required  functions  r,  s,  h^  and  h^  exist. 


When  these  cannot  be  found,  an  alternate  approach  is  that  of 
Fox  [9]  or  Lin  [13].  Again,  assuming  that  coefficients  of  (2.6)  are 
functions  of  u,  v,  e  only,  one  seeks  solutions  of  the  form 


u(£»n)  =  l  ek  u<k> (£,n) 

k=o 

(2.17) 

°° 

v(£,n)  =  l  ek  v  k  (£,n) 
k=o 


and  exnress 


(fl.S)  ■ 
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CO 

m  =  l  £k  (u,v)  , 

k=o 

(2.18) 


V  k  (k)  ,  . 

I  £  n±  (u , v ) 


00 

f(S)  =  l  ek  f (k) (?)  , 

k=o 

(2.19) 

g«)  -  I  ckg(k)(5)  • 

k=o 


Then  one  substitute  in  (2.6)  and  attempts  a  solution  to  the  successive 
problems , 


I  [m^X)  u<;j)(On)  +  n^x)  v^j)(C,n)]  =  0 
i+j  =k  n  n 


l  [m^l}  u?(j)(5,n)  +  nix)  vr(j)(C,n)]  =  o 


i+j=k 


2  £ 


with  initial  conditions 


( k)  s  _  -  \  r  ( k)  ,  w  V 

u  (5,0  =  f  (5)  , 


v(k)«,5>  =  g(k)(5)  , 


(2.20) 


(2.21) 


for  k  =  0,1,2, 


•  •  • 
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Once  we  have  u(£,p)  ,  v (C , n )  ,  we  can  attempt  a  solution 
of  (2.5).  We  substitute  these  expressions  in  P  and  Q  to  obtain 
PCCjPje)  >  Q(£>q>e)  which  may  be  in  series  form.  We  integrate  the 

equations 


x  =  P(£,n,e)t 

n  n 

=  Q(£,n>e)t^ 


with  boundary  conditions 


x(£,£)  =  £ 

t(£,£)  =  0 


(2.22) 


(2.23) 


from  the  initial  line  p  =  £  ,  along  the  curves  £  =  const.  and 
q  =  const.  respectively,  to  the  point  (£,n)  •  Thus 


n 

x(£,n)  -  £  =  /  P(£,a,e)  t^U.a)  da 


£ 

x(£,n)  -  n  =  /  Q(a,n,e)  t  (a,n)  da 

n 


(2.24) 


If  we  assume  that  we  let  e  0,  then  P  and  Q  are  constant;  we  have 
straight  line  characteristics  and  the  solution  is  known  exactly.  This 
suggests  that  we  attempt  a  power  series  in  e  as  did  Fox  and  Lin.  We 
express 


P(£,n,e)  =  l  ekP(k)(£,n) 

k=o 

(2.25) 

00 

Q(£,n,e)  -  I  ekQ(k)(C,n)  , 

k=o 


p  oxoesjcrxs  o«  itf  >  "is 


v 
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and  seek  solutions  of  the  form: 


x(?,n)  =  I  £k  x^  (£,n)  , 

k=o 

(2.26) 

00 

t(C,n)  =  l  ekt(k)a,n)  . 

k=o 


Initial  conditions  are  expressed  as 


x°U,S)  =  C 


7? 

'w' 

N 

vfn 

II 

o 

k  =  1,2,... 

(2.27) 

(k)(C,S)  =  o 

• 

• 

• 

• 

CM 

r\ 

1 — 1 

r> 

o 

ll 

We  substitute  in  (2.4)  and  assuming  that  the  series  converge  uniformly 


in  the  region  we  consider,  we  find 


00 

l 

k  (k)  .  v 

£  x  (5,n) 

00 

-5-1 

00 

l 

ek+J 

k=o 

k=o 

j=o 

CO 

l 

k  (k),  s 

e  x  (£,n) 

00 

-  n  =  l 

00 

l 

Ek+j 

k=o 

k=o 

j=o 

/  P(j)(£,<j)t^k)(£,a)da 

5 


(2.28) 


,(j) 


.00 

'a 


We  equate  coefficients  of 


k 

e 


and  assume  that 


and  Q° 


are  constants. 


Thus, 


x°(£,n)  -  £  =  P°t°(£,n) 


x°(£,n)  -  n  =  Q°t°(£,n) 


(2.29) 


I  -  * 


0  =  (2.3)a 


aw  ti9bianoo  sw  noigsa  9rf3  n± 
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n  k 


x(k)a,n)  =  P°t(k)(£,n)  +  /  l  P(n)(C,a)t(k  n)(£,a)da  , 

£  n=l  a 


(2.30) 


x 


(k) 


n  k 


(£,n)  =  Q°t(k)(£,q)  -  /  l  Q(n)(a,n)t^k  n)(a,n)da  . 


£  n=l 


By  addition  and  subtraction  we  obtain 


x°(£  ,n) 


n+i 

2 


+ 


,o 


P  +0 


n-£ 

2 


t°(£ ,n) 


(2.31) 


k 


(k) 


(£.n)  =  [Q0-P°]_1  /  ni± 


l  [P(n)(C,a)t^k  n)(C,a)+0(n)(a,n)t^k  n)(a,p)]da  C 


(k) 

for  k  =  1,2,...,  and  x  can  be  calculated  from  (2.30). 


In  the  problems  we  deal  with,  the  assumptions  made  on  P  and 
Q  are  satisfied  and  we  can  choose  an  upper  bound  for  q-£  ,  so  that, 
as  is  shown  in  Appendix  III,  these  series  for  x(£,q)  and  t(£,q) 
converge  for  all  required  values  of  (£,p)  ,  with  the  condition  that 

|A  +  B  +  A'  +  B*|  <_  |P°  -  Q°  ]  ,  (2.33) 

where  we  assume  that  for  the  values  of  £,£>0  considered:  jpj  <_  A  , 

|  Q  |  _<  B  ,  |P  |  <_  A’  ,  (Q  |  j<  Br  .  Such  conditions  involve  e  implicity. 
When  more  specific  information  is  available,  specific  bounds  on  e  and 
on  other  parameters  may  be  found;  for  example,  see  Fox  [9]. 

We  are  particularly  interested  in  solutions  at  values  of  t  where 


.32) 


lot  bnuod  i9qcai  ns  saoorio  nso  aw  bms  ba-Hellse  si 4  P 
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breakdown  occurs.  We  see  from  (2.31)  that  t  depends  mainly  on  n-£  , 
for  which  a  large  upper  bound  can  be  chosen. 


Once  we  have  the  solutions 


u  =  u(£,n)  , 
v  =  v(c,n;  , 
x  =  x(£,n)  , 
t  =  t(£,n)  , 


(2.34) 


then  with  an  assumption  that 
solve  for  £,n  in  terms  of 


3(x,t) 

9(C,n) 

x,  t  . 


1  0 

Hence 


it  follows  that  we  may 


u(x,t)  =  u(£(x,t),  n (x j t) ) 
v(x,t)  =  v(£(x,t),  n(x,t)) 


(2.35) 


expresses  solutions  for  Y  ,  Y  ,  from  which  Y  can  be  calculated.  We 

X  L 

deal  with  this  in  the  individual  cases. 


Kruskal 


A  slightly  different  approach  is  that  used  by  Zabusky  and 
[17].  They  assume  that 

P  =  P°  +  eP 


(2.36) 


Q  =  Q°  +  eQ 

where  P°  and  Q°  are  constants.  From  (2.22)  and  (2.23)  we  have: 

n  _ 

x(£,n)  -  K  =  P  t(£,n)  +e  /  P(£,a,e)t  (£,or)da 

K 

o  ,n  - 

x(£,n)  -  n  =  Q  t(5,n)  -e  J  Q(a,n , e)t  (a,n)do  . 


(2.37) 


•  BJm  rJenoo  9is 
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Substracting  we  obtain: 

-3  n  _  _ 

t(£»n)  =  [P°-Q°]  {(n-e)-e  /  [P(£,or,e)t  (C,a)+Q(a,n,e)t  (a,n)]da> 

'  o  o 


An  initial  approximation 


(2.38) 


(o) 


n  -  g 

P°  -  Q° 


(2.39) 


(k-1)  k 

is  assumed  and  substituting  t^  in  the  right  side  we  calculate  t 

Successive  approximations  for  k  =  1,2,...  are  calculated.  Convergence 

of  this  sequence  is  discussed  in  Appendix  II. 


If  instead  of  choosing  the  value  of  £  and  p  according  to 
the  characteristic  lines  x  intercept,  we  let  them  equal  their 
corresponding  Riemann  invariant,  £  =  r(£)  and  n  =  s(p)  ,  then  we 
can  write  (2.5)  as 


x  -  Pt  =0  , 

r  r 


x  -  Qt  =0 
s  s 


(2.40) 


Differentiating  with  respect  to  s  and  r  and  subtracting  we  obtain 


t  + 


P  t 
s  r 


rs  P  -  Q 


Q  t 

=  o 

P  -  0 


(2.41) 


an  equation  amenable  to  treatment  by  Riemann  function  methods  (see 
Zabusky  [15]). 


'  ■  - 


©dl)  /  olio  a  t  £.-9  vcf  .£  rr:  •  ;>3  s  c  •.  i-.n.' 


CHAPTER  III 


BREAKDOWN 


In  our  calculations  up  to  this  point,  we  have  discussed  solutions 


u,v  on  the  understanding  that  u  ,  v  ,  u  and  v  exist  and  are  continuous 

x  x  t  t 

in  the  required  region.  It  was  also  assumed  that  the  Jacobian 

J  =  -t  4-  0  there.  We  have  equations  in  the  form: 

3(C,n) 


x  =  P  t  x(C,0  =  K 

n  n 

X^  =  Q  t  t(S,£)  =  0 


(3.1) 


where  P  =  P(?,q,e)  and  Q  =  Q(£,q,e)  are  distinct  and  bounded  for  all 
values  of  £,  q,  and  we  assume  we  have  found  convergent  series  solutions 
for  (3.1),  in  the  form 


x(5  ,q)  =  l  £k  x^  (S  ,n)  > 

k=o 

(3.2) 

t(c,n)  =  l  e  t  (S,n)  • 

k=o 


We  can  now  express  the  Jacobian: 


J  =  x_  t  -  x  tr 

S  n  n  £ 


as 


J  =  [P  -  Q]t?  t^ 


(3.3) 


' 


o=>I 
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If  it  should  happen  that  J  =  0  at  some  point  (£,ri)  then 
clearly  our  transformation  from  the  (x,t)  to  the  (£,ri)  plane  is 
singular  and  either 

t_  =  x„  =  0  ,  or  t  =  x  =  0  , 

£  £  n  n 

in  which  case  £  ,  £  or  lr  >  h  are  unbounded.  Thus  u  =  u„£  +  u  p  , 

t  x  t  x  x  £  x  r)  x 

u  .  v  and  v  become  unbounded  and  we  cannot  admit  u,v  as  solution 
t  x  t 

to  our  problem  at  such  points.  Points  or  curves  along  which  such  a  break¬ 
down  occurs  define  a  boundary  of  the  domain  in  which  we  have  found  a 
solution  to  our  initial  value  problem.  In  this  Chapter  we  examine  the 
occurrence  of  such  curves  and  the  resultant  "folding"  of  the  characteristic 
plane . 


From  (3.3),  we  see  that  a  locus  of  points  in  the  (£,n)  plane 
where  J  =  0  is  obtained  from  the  curves 


t  (£>n)  =  0  , 

or 

t  (£,n)  =  0 
n 


If  we  approximate 


t,  =  I  ek  £'-(k)(?,') 

^  k=o  ^ 


0  , 


(3.4) 


(3.5) 


by  taking 

tj?(£,n)  +  e  t5(1)(£,n)  =  0  ,  (3.6) 

then  we  will  approximate  points  at  which  (3.4)  is  satisfied.  The  assumption 


here  is  that  we  have 


hy, 
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00 


|  l  ekt^(C,n) |  <  K  e 
k=  2  ^ 


(3.7) 


for  some  constant  K  ,  in  the  neighbourhood  of  points  which  satisfy  (3.6). 
Thus,  if  e  is  small,  the  solution  to  (3.6)  cannot  differ  significantly 


from  that  of  t^(£,n)  =  0  .  The  analysis  is  similar  for  t^(£,n)  =  0  . 


We  have,  as  yet,  been  unable  to  form  a  general  theorem  under 


which  our  assumption  (3.7)  can  be  shown  to  hold.  It  would  appear  that 
we  must  consider  each  problem  individually.  If  one  refers  to  Zabusky, 
where  a  specific  Riemann  function  is  known,  and  hence  a  specific 
representation  of  the  solution  t(£,n)  it  is  possible  to  verify  these 
assumptions.  One  such  verification  is  considered  in  Appendix  4. 

As  will  be  shown,  in  the  problems  we  consider,  it  is  because 
we  can  write 


P(1)(£,n)  =  h1(C)  +  gl(C,n) 


(3.8) 


Q(1)  (S,n)  =  h2(n)  +  g2(s,n) 


where  h^,  g^  are  bounded,  and  h^  ,  i  =  1,2,  are  non1- const  ant , 
that  singularities  occur.  To  show  this  we  assume  that  P°  and  Q° 


are  distinct  constants  and  that 


00 


k=o 


k=o 


From  (2.32)  we  see 


83n±og  'io  boorfit'oJj';; 

j  r.'XC  j.  o  .  /"i 


V 
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=  [Q°-P°]  1|^/ri  [P(1)(^a)t°(c,a)+q(1)(a,n)t°(a,n)]da 


=  [Q°-P°]  2  {-  /  Pru;(?,a)dc  +  Qu;a,n)-PU;(?,n)) 


(1) 


(1),,  „(1) 


[P°-Q°]  2  (-  /n  |y  h1(C)da  -  /n  gia,o)da  +  P(1)(^,n)-Q(1)  (C » n ) > 


[P°-Q°]  2  (-( n,S)h[(£)  +  B(c,n)> 


(3.10) 


where  B(£,q)  represents  terms  uniformly  bounded  in  £,q  ,  h|(£) 

represents  -r—  h..  (?)  and  we  have  assumed  that  f  '  2  (£,a)dcr  is  uniformly 

d?  1  ?  d?  61 

bounded  in  £,p  . 


Thus,  when  we  approximate  breakdown  by  solving  the  equation 


0  =  t°  +  e 


(3.11) 


we  have 


0  =  [P°-Q°]  2  {Q°-P°-e(n-5)h'(5)}  + 


or 


0  =  P°-Q°  +  e(n-0  h|(?)  4-  0(e) 


(3.12) 


In  view  of  (3.7)  we  are  in  fact  setting  t^  =  0  4-  0(e)  so  that,  in 
setting 


n 


=  n(S)  "  S  ~'lhr(Ti  ’  e  e  {o|(P0-Q°)h^(o)  <  0}  ,  (3.13) 


0 


. 
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we  are  approximating  the  curve  at  which  t  =  0  .  A  different  approach 
to  the  study  of  the  singularities  is  given  by  Zabusky  and  Kruskal  [17] . 
Since 


3  (x ,  t)  _  3(x-P°t>  x-Q°t)  #  _ 1_ 


8(5, n) 


8(5, n) 


[P°-Q°] 


(3.14) 


where  P°  and  Q°  are  distinct  constants  and 


P  =  P°  +  eP  ,  Q  =  Q°  +  eQ  , 


we  can  consider  a  different  Jacobian.  The  derivatives  involved  here 
express  the  rate  of  deviation  from  the  linear  characteristics  on  which 
we  are  perturbing.  If  we  write  (2.36)  in  the  form 


x  -  P°t  =  £  +  ® 
x  -  Q°t  =  n  +  © 


(3.15) 


where 


and 


n  _ 

®  =  I  P(5,o,e)  t^(£,a)  da 


n 

©  "  /  Q(5,n,e)  t^(a,n)  da 


(3.16) 


then  the  new  Jacobian  has  the  form 


a  +  ®t)d  +  ®n)  -  ®n©5 


(3.17) 


©n 


is  of  order 


2 

£ 


an  approximation  to  the 


Since 


®  +  n  3 
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location  of  zeros  of  J  should  be  provided  by  setting 


1  +  ©  =  0  +  0(e2) 
1  +  @  =  0  +  0(e2) 


(3.18) 


The  procedure  is  now  essentially  the  same  as  above. 

Finally,  a  prediction  of  the  time  when  singularities  occur  can 
be  arrived  at  without  actually  solving  the  problems.  The  method,  which 
involves  examining  the  derivatives  r  ,  s  of  the  Riemann  invariants,  was 

X  X 

applied  by  P.  Lax  [12]  to  the  problem 


Y  =  — r  (1  +  Y  )  Y 
tt  2  x  xx 

c 


Y(x,0)  =  b  sin 


t  >  0 


(3.19) 


to  predict  a  breakdown  time,  t^  ,  corresponding  to  those  obtained  by 
Zabusky  and  Kruskal. 


Consider  a  set  of  Riemann  invariants  r  and  s  of 


Y  =  K  Y 
tt  xx 


(3.20) 


and  let  f '  and  fN  denote  differentiation  of  a  function  f  along  the 
characteristic  directions.  Let  P,Q  be  distinct  characteristic  values  of 
the  matrix  A  in  (2.1).  Choose  h  such  that: 


h  =  P  / P-Q 
s  s 


p  -i£ 

s  8s 


(3.21) 


If  we  set 


z  =  r  e 
x 


(3.22) 


Y(  y  +  £)  r- 


~  HX2  d 


2nifanoqa9'iiOD  « 


i  aa.xfiiievnx  an  ni9i.fl  do  39 3  &  :c9bianoD 
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a  -  -e  h  Pr  ,  (3.23) 

gp 

where  P  =  —  is  assumed  negative,  we  obtain 
r  dr 

7/  +  a  7}  =  0 

Two  theorems  are  given  in  Lax’s  paper  [12]: 


"Theorem  1.  Let  z(t)  be  the  solution  of  the  initial  value  problem 

-^7  =  a(t)z2  z (0)  =  m  (3.24) 

in  the  interval  [0,T]  .  Suppose  that  the  function  a(t)  satisfies  the 
inequality 


and  suppose  that 


0  <  B  <  a(t)  0  _<  t  <_  T  , 

m  is  positive,  then  T  <  (mB) 


Theorem  2.  Suppose  that  a(t)  satisfies  the  inequality  J a ( t ) J  <  C  ; 
then  the  initial  value  problem  (3.24)  has  a  solution  for  ] t  j  <  |mC|  ^  ." 

If  we  assume  that  there  is  a  positive  maximum  value  ,  m^  ,  of 
z  on  the  initial  line,  and  provided  the  bounds  B  and  C  exist,  then 
the  theorems  to  determine  bounds  for  the  time  t^  ,  namely, 

(Cm^)  1  <  t  <  (Bm^)  ,  (3.25) 

at  which  z  is  unbounded.  Because  z  =  e  r  ,  and  h  is  bounded,  r 

x  x 

is  also  unbounded  at  that  time. 


sxtt  ,  n  vj.3  1  a  q 


vrillsupani  »rft  aai  alt  .6  (3)b  36,1.  aaocp  v? 

3  e  art  ,  e,r  rod  9  oj  ; 
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To  see  how  this  applies  to  the  present  problem  we  let  K  =  —  (1+eY  ) 

c 

in  equation  (3.20)  and  as  in  Chapter  II  we  write  the  system  in  terms  of 

u  =  Y  ,  v  =  Y  .  The  characteristic  values  P,  Q  are  -K,  K  and  we 
x  t 

find  Riemann  invariants  r,  s  such  that  r  =  K,  r  =1,  s  =  -K,  s  =1, 

u  v  u  r 

so  that 


r  =  v  +  L(u)  ,  s  =  v  -  L(u)  ,  (3.26) 


where  L(u)  =  K  .  In  order  to  compute  P^  ,  we  differentiate 

r  -  s  =  2  L (v)  ,  (3. 27) 


to  obtain 


1=2  L(u)  v  =  2  K  u  (3.28) 

u  r  r 


Thus 


P  =  P  u 
r  ur 

-K 

u 

2K 

1 

=  ~  (1+  Y  )  2  ,  |Y  |  <  -  ,  (3.29) 

4  x  ’  1  x1  e  * 


Further , 


r  =  v  +  L  u  =  v  +  K  u 
x  x  u  x  x  x 


so  that  the  initial  values  are 


1 

r  (0)  =  -  (1+eY  (0))2  Y  (0) 

X  C  X  XX 


(3.30) 


Here  (0)  denotes  the  initial  point  of  a  given  characteristic  line. 
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Under  the  assumption  that  initial  conditions  have  small  variation,  we  can 
assume  that  the  quantities  and  h  vary  little  from  their  initial 

values  P  (0)  and  h(0)  ,  so  that  one  can  estimate: 


(Cm.)  ^  =  [max  (-e  p  (o)),  max  (r  (0)  e*1^)]  ^ 

j-  r  x 


[max  (-P  (0))«max  (r  (0) ) ] 

L  X 


-1 


i  (l+eM)2 
e  n  bn. 2  ^  r 

4  U  '  c 


n-1 


bn' 


£2J 


(3.31) 


If  we  assume  that  e  is  small,  we  have 


(Cm1> 


-1  4  i 


r 

en  be 


(3.32) 


Similarly  (Bm^)  ^  has  the  same  value,  so  that  we  have  a  prediction  of 

2 


breakdown  time  t,  = 

b 


4Jt 


cben' 


1  o 

But  for  K  =  —  (1+eY^) 
c  x 


p  =  2, 3, . . .  ,  we  have 


P  = 
r 


2  VP“1 

£  P  Yx 

4 (l+cYP 

X 


(3.33) 


=  0  for  Y  =  0 
x 


so  that  there  exists  no  lower  bound  B  for  |e  P  j  and  the  theory  does 
not  provide  an  upper  bound  for  the  interval  of  existence.  Thus  breakdown 
is  not  predicted.  For  a  more  detailed  investigation  of  the  correspondence 
between  the  (£,n)  and  the  (x,t)  planes  in  the  neighbourhood  of  the  singular 


' 
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lines  we  can  refer  to  theorems  due  to  Cragg  [5],  formulated  for  the  study 
of  the  breakdown  of  the  hodograph  transformation. 


From  [5]  we  shall  note  two  theorems.  Let  us  first  of  all  define 
a  singular  line  as  a  line  along  which  J  =  0  .  Denote  it  by  S  in  the 
(5>n)  plane  and  by  L  in  the  (x,t)  plane.  We  define  a  singular  point 
of  order  n  to  be  a  point  at  which  all  derivatives  of  J  of  order  less 
than  n  are  zero  as  well  as  J  .  We  discuss  the  singular  line  S  , 
corresponding  to  t  =  0  . 


We  quote  from  [5],  where  we  have  replaced  the  symbols  used  there 
by  the  corresponding  ones  as  used  in  this  thesis. 


’’Theorem  1.  Let  a  curve  S  in  the  (£,n)  plane,  each  point  of  which 
is  singular  of  order  one,  correspond  to  a  curve  L  in  the  (x,t)  plane. 
Then,  save  in  exceptional  cases,  with  each  point  of  S  is  associated  a 
privileged  direction,  such  that  to  any  curve  crossing  S  in  that  direction 
there  corresponds  a  curve  having  a  singular  point  (in  the  usual  geometrical 
sense)  on  L  ,  while  to  any  curve  crossing  S  in  a  direction  other  than 
the  privileged  direction  there  corresponds  a  curve  touching  L  . 

Thoerem  2.  If  a  cusp  of  a  curve  is  defined  as  a  singularity  such  that 
the  adjacent  parts  of  the  curve  on  either  side  of  it  have  the  same  tangent 
but  in  opposite  senses,  then  (i)  except  where  the  singular  line  S  is 
itself  in  the  privileged  direction,  all  curves  crossing  S  in  the  privileged 
direction  correspond  to  the  curves  cusped  on  L  ,  and  (ii)  when  the 
singular  line  S  is  itself  in  the  privileged  direction  at  a  point  and, 


",  !  cyj  t  " 
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in  the  notation  of  Theorem  1, 
is  cusped  ... 


t  +  f" (5)t  4  0 

KK  a 


the  limit  line  L 


Theorem  3.  All  points  in  the  vicinity  of  a  singular  line  S  of  order 
one  correspond  to  points  on  one  side  of  the  limit  line  L  . M 


In  the  development  of  these  theorems,  it  is  shown  that, 
corresponding  to  curves  where  t  =  x  =  0  ,  the  curves  p  =  const, 
are  not  in  general,  in  the  "privileged  direction"  and  thus  tangent  to  L 
in  the  (x,t)  plane,  also  that  L  is  cusped  in  the  (x,t)  plane  at  points 


where  t  =  0  .  The  parameter  a  ,  generally  different  from  p  or  £» 

3  x 

that  for  which  —  =0  at  a  given  singular  point  under  consideration,  and 

da 

f  is  defined  such  that  a  =  f(£)  is  the  singular  curve  S  . 


From  (3.12)  we  see  that  t  =  0  +  0(e)  when 

A 

e(p-C)  h  (£)  =  0  .  (3.34) 

d£ 

which  corresponds  to  the  minimum  of  p(£)  ,  (3.13),  where  —  p(^)  =  0  . 

d£ 

We  do  not  consider  points  where  singular  lines  interesect  in  the  (£,p) 
plane,  that  is,  where  t_  =  t  =  x  =  x  =0. 

F  £  p  Z  p 

Consider  Figures  V  and  VI.  Theorems  1  and  2  imply  then  that  the 
initial  breakdown  point  G  ,  the  extremum  of  the  t^  =  0  line,  where 
t^  =  0  ,  maps  into  a  cusp  in  the  (x,t)  plane.  The  lines,  p  =  const., 
and  E,  =  const,  are  tangent  to  L  at  points  such  as  A,  B  and  F  . 
Theorem  3  implies  that  points  C  and  D  should  be  on  the  same  side  of 
L  ,  and  similarly  for  D  and  E  . 


is 
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G 


LINES  CROSSING  S  IN  THE  (5, n)  PLANE 

Figure  V 


IMAGE  LINES  CROSSING  L  IN  THE  (x,t)  PLANE 


Figure  VI 


' 
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From  this  comes  the  notion  of  the  "folding"  of  the  characteristic 

plane. 


FOLDING  OF  (£,n)  PLANE  ABOVE  CORRESPONDING  POINTS 

OF  (x,t)  PLANE 

Figure  VII 

The  folded  sheet  (Figure  VII),  represents  the  (£,n)  plane  placed  above 
the  corresponding  points  in  the  (x,t)  plane.  For  each  point  of  the  (£,n) 
plane  there  is  a  unique  solution  (u,v)  .  Thus  where  the  "fold"  occurs 
there  are  three  values  of  (£,n)  corresponding  to  one  point  (x^,t^)  . 

This  leads  to  no  contradiction  of  the  existence  theorem  [4],  as 

only  one  of  those  three  points  lies  between  the  singular  line  and  the 

initial  line  in  the  (£,n)  plane.  While  there  are  three  values  of  (u,v) 

for  (x^jtj,)  inside  the  region  of  existence  of  a  solution,  only  the  value 

corresponding  to  the  point  in  the  (£,n)  plane  between  the  initial  line 

and  the  singular  line  can  be  considered  a  solution.  Those  points  (£,n) 

in  the  "shadow"  of  the  singular  line  are  not  in  the  domain  of  existence 

of  the  solution  because,  contrary  to  the  existence  theorem,  u  ,v  ,u  ,v 

X  X  U  L 
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are  unbounded  in  the  domain  of  dependence  of  the  point.  Thus  the  upper 
boundary  of  the  domain  of  existence  consists  of  singular  lines,  (e.g.  EB 
in  Figures  VIII  and  IX)  and  characteristic  lines  running  up  from  the  vertex 
of  a  cusp  (e.g.  BF)  . 

It  will  be  seen,  when  we  consider  the  particular  problems,  that 
the  conditions  t_v  =  0  and  t  =0  hold  at  the  initial  breakdown 

KK  nn 

points  and  that  the  shape  of  t^  =  0  is  similar  to  that  in  Figures  V  and 
VII.  The  characteristics  of  Figures  VI,  VII,  and  IX  are  exagerated, 
compared  to  the  numerical  results  performed  for  the  fixed  end  problem. 


EXAMPLE  OF  DOMAIN  OF  EXISTENCE  IN  (£,n)  PLANE 


Figure  VIII 


* 
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EXAMPLE  OF  CORRESPONDING  DOMAIN  OF 
EXISTENCE  IN  (x,t)  PLANE 


Figure  IX 


CHAPTER  IV 


THE  FIXED  STRING  PROBLEM 


In  this  chapter,  we  consider  a  generalized  fixed  string  problem 


Y  =  (1  +  eYP)q  Y 
tt  x  xx 


0  <  x  <  1 

t  >  0  , 


(4.1) 


subject  to  the  boundary  conditions, 


Y(x ,  0)  =  a  sin  nllx  , 


Yfc(x,0)  =  0 


;  0  <  x  <  1  , 


(4.2) 


Y(0 , t)  =  Y(1 , t)  =  0  ,  t  >  0  , 


where  n  and  p  are  positive  integers  and  q  is  a  non-zero  integer. 

We  will  of  course  require  that  1  +  eYP  >0  in  order  to  ensure  boundednes 
in  the  various  cases  we  consider. 


We  write  the  system  as 


u  \ 

0 

-1 

f  u 

l°\ 

+ 

1  v  1 

i-$2 

0 

,  V 

“U 

x 


t  >  0 


-OO  <  X  < 


(4.3) 


with  initial  conditions 

u  =  U(x,0)  =  anil  cos  nllx  ,  -oo  <  x  <  °° 
v  =  V(x,0)  =  0  , 


(4.4) 


where 
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<Ku)  =  (1  +  eup)q/2  . 


(4.5) 


It  will  be  shown  later  that  this  periodic  extension  of  initial  data 
satisfies  end  conditions.  We  bound  a  and  e  so  that  0  <  c(nan)P  <  1 
this  will  mean  that  (1  +  euP)  >  0  ,  which  provides  for  the  boundedness 
of  $(u)  and  its  derivatives. 

Eigenvalues  for  the  matrix  are  P  =  -$  ,  Q  =  $  ,  with  left 
eigenvectors  (m^,n^)  =  ($,1)  and  (m  j^)  =  (-$,1)  .  As  in  Chapter  II, 
we  express  the  system  in  terms  of  characteristic  variables.  Thus  we  have 


-$u^  +  v^  =  0 

$u  +  v  =  0 

0  'n 


x^  +  $t  =  0 
x  -  $t  =  0 

n  n 


5 


with  boundary  conditions, 


u(s,s)  =  f(s)  , 
v(£,€)  =  0  , 


xtt,a  =  K 

t(£,0  =  o  , 


where  f(£)  =  anil  cos  nIT£  • 


(4.6) 


(4.7) 


(4.8) 


We  employ  the  Riemann  invariants  to  solve  for  u(£,n)  ,  v(£,p) 

We  select  hp  =  "  ^  ’  h2  =  2^  S°  that 


u 


$(u) 


s  -  —  <Ku) 
u  2 


1 

r  =  -r 
v  2 


s  =  - 
v 


(4.9) 


1 

2 


-  41  - 


The  Riemann  invariants  are: 


u 


r  =  r(?)  = 


4  /  $(a)da  +  ^ 


o 

u 


s  =  s ( n )  =  *9  /  <KcOda  - 


v 


(4.10) 


Let  us  consider  first  the  simpler  case  where  p  =  1  ,  q  ^  -2  ,  then, 


u 


r  +  s  = 


J  $(a)da 
o 


[e(f  +  1)] 


-1 


£ 

[(1  +  cu)2  -  1] 


(4.11) 


r  -  s  =  v 


In  general  if  we  assume  that  |u|  <  -qy  then  we  have 

^ r ’  S\  =  4$^  0  and  we  can  invert  the  expressions  so  that,  in  this  case 
3(u,v)  2 

with  p  =  1  , 


In  1/(1+  'o’) 

u  =  u(r+s)  =  —  [e(-+  +  1)  (r+s)  +  1] 

£  Z 


1 

£ 


(4.12) 


v  =  s  -  r 


Under  the  assumption  that  we  are  dealing  with  a  region  where 
u , v  exist  and  are  properly  bounded,  we  have^*r  =  r(£)  and  s  =  s(q),  constant 
along  characteristic  lines,  and  whose  values  we  obtain  from  the  initial 
line  : 


s(a)  =  r (a)  = 


[2c  (|  +  1)] 


1  +  3- 

1  •  [(1  +  cf (a) )  2  -  1] 


(4.13) 
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Combining  (4.12)  and  (4.13)  we  obtain: 

u(S,n)  =  u(r  (n) ,  s(O)  , 
v(?,n)  =  r(0  -  s(n) 


(4.14) 


have 


Because,  as  we  see  in  (4.9), 


r  and  s  are  positive,  we 
u  u 


max  r  =  r(u^) 

(4.15) 

max  s  =  s(u^) 


where  =  max  u  ,  thus 


max  (r  +  s)  =  r(u^)  +  s(u^) 


(4.16) 


If  we  set  r+s  =  y  and  consider  (4.11),  we  see  that  4^-  =  $(u)  >  0  and 

dy 

thus  =  $(u)  ^  >  0  ,  for  properly  bounded  values  of  u  .  Because  u 
is  monotonically  increasing  in  r  +  s  ,  max  u  =  u(max  (r+s))  = 
u(r(u^)  +  s(u^))  .  Thus  the  maximum  value  of  u  corresponds  to  the 
maximum  value  u^  along  the  initial  line  and  thus  u  <  ^  By  similar 

analysis  we  have  a  lower  bound.  Our  solution  u(£,p)  in  (4.14) 


satisfies  the  condition  u  <  —  for  all  values  of  £  and  p 

£ 


We  proceed  with  the  general  case  where  we  will,  in  addition, 
have  to  resort  to  series  expansions.  Consider 


■ 


7  XW: 


43 


u 


r+s  =  I  (1  +  euF) 
o 


Pxq/ 2 


u 


=  /  {1  +  e  aP  +  ^  (^  -  1) 


Pn2 


(coO 

2! 


+  ... }da 


(4.17) 


p+1  2  2p+l 

=  u  +  5.  _  +  £(£_!)  £_.-P. _ _  + 

2  p+1  2  ^  ;  2!  2p+l 


/  2 

Because  of  our  bound  eju^]  <  1  ,  (1+eu^)^  is  an  analytic 

function  of  u  ,  and  so  is  its  integral.  Thus  the  series  (4.17)  converges 
and  has  a  non- zero  derivative  at  u  =  0  .  As  a  result,  it  can  be  inverted 
to  give 


u 


=  l  ai(r+s)1  , 


1=0 


for  some  bound  6,  r+s  <  6  .  The  first  few  terms  are: 


u  = 


(r+s)  ’  Icf+oJ  (r+s)P+1 


+  £ 2  r — 3 - q(q-2>  ]  (r+s) 2p+1  + 

L4(p+1)  8(2p+l)U  s;  * 


(4.18) 


In  order  to  determine  the  region  of  convergence,  we  let 


y  =  r+s  and  consider 


u 


y  =  /  (1  +  ea^)^^  da 

o 


(4.19) 


so  that 


du 

dy 


,  p. q/2 

(1+eu  ) 


u(0)  =  0 


(4.20) 
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Now,  consider  the  differential  equation 


du 

dq 


=  f(u) 


(4.21) 


u(0)  =  0 


where  f(u)  is  analytic  for  |u|  <  b  and  |f(u)|  <  M  there.  It  is  well 
known,  for  example  see  Goursat  [11],  that  the  series  solution  for  this 
equation  is  dominated  by  the  series  solution  for 


dU 

dy 


M 


U(0)  =  0 


Further,  the  solution  of  (4.19)  is  given  by 


u  -  jr  =  M 
2b  y 


(4.22) 


U  -  2bU  +  2bM  =  0 

y 


U  =  b [1  + 


C 


2M 


(4.23) 


which  is  analytic  in  y  for  | y |  <  —  .  In  our  case 

_  SL 

M  =  max  (1  +  eu^)  2 


= 


_  £ 

/'i  b  \  2 

(1  -  EUp 

_  4 

(1  +  euP)  2 


q  >  0 

q  <  0 


(4.24) 


where  u^  =  max  u  ,  and  b  <_ —  .  Thus  (4.18)  is  convergent  for 


- 1) 

»S»vnoo  *1  (81. A)  auifT  .  te  ,  „  ZBln  .  „  M(Mh, 
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r+s  < 


_  £ 

(1  -  £UP)  2 

to? 


q  >  o  , 


r+s  < 


(1  +  euP)  2 

2? 


(4.25) 


q  <  0 


Clearly,  we  may  choose  e  small  enough  to  cover  all  values  of  r+s  . 


Again,  as  r  and  s  are  constant  along  lines  £  =  const.  and 
n  =  const.,  respectively,  they  have  values  from  the  initial  line  calulated 
as  in 


r(a) 


s  (a) 


£ 

(1  +  eaP)2 


da 


[ f (a)  +  £ 


a  f(a)P+1 

2  p+1 


+  .  .  .  ] 


(4.26) 


Again  r  and  s  ,  and  thus  u(y)  =  u(r+s)  reach  their  maximum 
and  minimum  values  where  the  initial  function  u  =  f(a)  reaches  its 
maximum  and  minimum.  Thus,  substituting  the  results  from  (4.26)  into 
(4.12)  and  (4.18),  that  is,  the  values  for  r($)  ,  s(n)  into  v  =  s-r 
and  u  =  u(r+s)  ,  we  obtain  u(£,n)  ,  v(£,n)  and  |u(£,ri)  I  <  max  f  (a)  <  — 

_  p 


We  next  use  the  method  of  Fox  and  Lin,  outlined  in  Chapter  II, 
to  obtain  first  terms  of  the  expansion  of  the  solution  for  x  and  t  . 
Since  we  have  | euP j  <  1  ,  we  can  express 


$  =  (1  +  £uP ) 2 


=  1  +  -|  £UP  + 


£  (£  ~  1) 
2  k2  j 

21 


(euP) 2  + 


•  •  • 


(4.27) 


- 
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and  from  (4.18)  we  have 


u  =  (r+s)  - 


2 (p+1) 


(r+s)^+^  +  .  .  . 


so  that 


$  =  1  +  e{  (r+s)P  -  2(p+l)'(-r+s-)2P  +  *  •  *  ^ 


+  ^21  £2  {(r+g)2p  +  _  } 


1  +  e  (r+s)P  +  e2 [-ill— (r+s)2P+  .  .  . 


Since 


f  (a) 


(a)  =  s(a)  -  ^  J  (1  +  eaP)P//2  da 


.  1  f /  n  ,  eg  f(a)P+1 
2  f  a  4 (p+1) 


where  f  (x)  =  anil  cos  nllx  ,  we  have 


r+s  = 


|(f(n)  +  f(O)  +  4(p+1y(f(n)P+1+f  (Op+1)+. 


•  • 


finally , 


»(C,n)  -  1  +  ~rrr  (f(Op+1+f(n)p+1) 

0P+1 


,  2  { fq(q-2)  _  .  -.q^P-.]  (f(g)  +  f(nK 2p  + 

1  L  8  4 (p+1)  ^  2  ; 


+  8(fe  (f<°  2  (f(5)P+1  +  £(^P+1» 


(4.28) 


(4.29) 


(4.30) 


(4.31) 


+  0 (e3) 


(4.32) 


. 
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We  now  set  P  =  $  and  Q  =  -$  so  that,  from  (2.31) 


of>.  \  n  +  S 
x  (S,q)  =  — 2 — 


t°(S,q) 


=  n  ~  A 
2 


(4.33) 


and  from  (2.32) 


t(1)  -  -  \lf  i-f(^i)±^>)Pda  +  f 


da] 


=  — ^  l  C  )(f(C)p  k  +  f(n)p  k)  ^  f(0)k  da  , 

2P  k=o  R  £ 


(4.34) 


x 


(1)  -  [  (P)(f(5)p-k  -  £(n)p-k)  /”  f(a)k  do 

2P  k=o  ^  E, 


In  particular,  for  p=l,  q=l,  n  =  1,  f(S)  =  all  cos  Ilf  we  have: 


Cl')  alt  a 

£  '(S,q)  =  Y^(rr?)(cos  JIE-  cos  fin)  +  -g(sin  fin  sin  fiS) 


x^(s,q)  =  —(n-C) (cos  ns  -  cos  nn)  > 


and  for  p  =  2 ,  q  =  1,  n=l  we  have 


(4.35) 


t<'1'>(s,n)  =  -  ^ ip-  (n-s) (cos2ns  +  cos2nq) 


a2n 

16 


(cos  US  +  cos  fin) (sin  Hq  -  sin  US) 


a2U  1 

■jj-  [  Cn— C)  + -^(sin  2nq  -  sin  2HS)] 


(4.36) 


x^(s,q)  =  ^!p .  (n-S)(cos2ns  -  cos2nn) 

a2n 

-* - (cos  ns  -  cos  Hn)  (sin  Hq  -  sin  nS) 


16 
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In  general,  when  we  set  f(x)  =  (an-H)  cos  nHx  ,  and  observe  that 


Gk(o)  -  /  cosk  ado 

;  nil 


1  J  1  k-1  .  k-1  1  k-3  ,, 

=  — r  <  7-  cos  a  sin  a  +  — —  - — -  cos  a  sin  a  (4.37) 
nil  I  k  k  k-2 


f  k-1  3/0.1  s  ] 

k  ’ ‘ *  4  z  +  2  *  Sin  °  COS  a'  k  even 

+  ...  +  < 

>  > 

k-l.  2 

■  .  .  .  .  ~  sin  a  k  odd 

k  3  J 

k 

/ 

+c 


Then  we  obtain: 


t(1)(£,n) 


-  ■^3  (ann)P  l  [(P)(cosP  Vn£  +  cosP  Vim) 
2P  k=o 


(Gk(nnn)  -  Gk(nIlO)]  • 


(4.38) 


x(1)(£,n) 


— p'+3  (anII)P  l  [(P)(cosP  Vll£ 

2?  k=o 


cosP  Vim) 


(Gk(nnn)  -  Gk(nn5))]  . 


(4.39) 


Taking  derivatives  of  these,  one  could  calculate  further  terms 
of  the  solution  series  according  to  (2.32). 

Under  our  assumption  that  we  deal  with  a  region  where  J  ^  0  , 
we  attempt  to  find  £(x,t)  ,  n(x,t)  •  We  seek  expressions  of  the  form 


£  =  £°(x,t)  +  e  £  ^  (x , t ) 


n  =  n°(x,t)  +  e  (x,t) 


(4.40) 
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The  terms  £°  >  9°  a re  the  unperturbed  linear  expressions  x  -  t  ,  x  +  t  . 

From  expressions  of  x  and  t  , 


t 


+  e  t 


(1) 


x 


n  +  g 

2 


+  £  X 


(1) 


(4.41) 


we  obtain 


£  =  x  -  t  +  e(x^^  -  t^^)  , 


x  =  x  +  t  4-  e  (x^'L')  +  t^) 


for  p=l,  q=l,  n  =  1  ,  from  (4.35),  we  have 


an  3 

£  =  x-t  +  e[—  (n-£)  cos  Up  +  -g  (cos  Tin  -  sin  H£)] 


ail  a 

=  x-t  +  e  [—  t  cos  H(x+t)  +  —  cos  Hx  sin  nt] 


(4.42) 


all  a 

n  =  x+t  -  £  [n —  t  cos  n  (x-t)  4-  —  sin  Hx  cos  lit] 

4  4 


For  t  small,  the  series  of  which  these  are  the  first  terms  surely  give  a 
good  approximation  where  £  is  small.  But  as  t  ->  t^  ,  higher  order  terms 
may  become  significant. 

Similarly  for  p=2,  q=l,  n  =  1  ,  from  (4.36),  we  have 
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cos  Il(x-t)  cos  Hx  sin  nt  +  7-  cos  2IIx  sin  2Ht 


(4.43) 


n  =  x  +  t  -  e  — t  cos2n(x+t) 


a  'll 

-  e  — r—  cos  n(x+t)  cos  nx  sin  nt 


Thus,  combining  values  of  u(r,s)  ,  v(r,s)  ,  values  of  r(£)  ,  s(n)> 

and  those  of  £(x,t)  ,  p(x,t)  ,  one  can  arrive  at  u,v  =  (U(x,t),  V(x,t)) 
form  which,  Y(x,t)  could  be  calculated. 

We  have  considered  the  problem  in  the  half  plane,  but  the 
solution  also  satisfies  the  end  conditions  of  the  fixed  end  problem: 


Y(0 , t)  =  Y(1 , t)  =  0 


(4.44) 


x 


1 


0 


1 


x 


1 


CHARACTERISTIC  LINES  INTERSECTING  IN  x  =  0  . 


Figure  X 


51  - 


Because  U(x,0)  =  U(-x,0)  and  V(x,0)  =  0  , 

from  (4.10)  we 

have  r (a)  =  s(a)  =  s(-a)  •  Thus,  from  (4.11) 


v(c,-a)  =  r(a)  -  s(-a) 

(4.45) 

=  0 

and 


u(a,-c)  =  u(r(a)  +  s(-a)) 

=  u(r(-a)  +  s  (a) ) 

(4.46) 

=  u(-a,o) 


and  further 

$(a,-a)  =  $(-a,a) 

(4.47) 

From  this  we  see  that  the  slopes,  P  =  $  and  0  =  -$ 

of  the  characteristic 

lines,  and  the  lines  themselves  are  symmetric  about  x 

=  0  .  It  then 

follows  that  for  any  point  x(£,p)  =  0  ,  we  have  £  = 

-p  and  so 

V(0,t)  =  v(-n,n) 

=  0 

(4.48) 

Thus,  integrating  along  x  =  0 


C1 

Y(0,t1)  =  /  V (0 , t)dt 

o 

(4.49) 

=  0 

A  similar  analysis  can  be  performed  for  x  =  1  as  initial  data  is  also 
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symmetric  about  x  =  1  . 


We  can  now  indicate  possible  breakdown  times,  using  the 
approximations  we  have  to  t  .  We  do  not,  in  fact,  prove  that  breakdown 
occurs,  but  that  only  under  the  assumption  that  the  neglected  terms  are 


small,  we  calculate  points  where  t  „  or  t  are  zero.  From  (4.34), 

S  n 

we  have  that 


t 


(1) 

K 


U  >n) 


p 

-9-  y 

2P+3  k=o 


P 

—9-  y 

2P+3  Jo 


(P)  [f(C)p  k  +  f (n)p  k]f(5)k 

(p)  (p-k)f (0P_k_1  f'(C)  /  f(o)kdo 


(4.50) 


and 


t(1)(C,n)  =  -t(1)(n,0 

n  4 


(4.51) 


We  observe  that 


is  0(1)  only  in  those  terms  for  which 


p  u 

f  f (a)  do  =  (n-C)h(C,n)  >  (4.52) 

5 

in  which  case  h  is  constant.  In  this  case,  e(n-£)  is  0(1)  when 
(n-5)  -  7  • 

Setting  f(x)  =  anil  cos  nllx  in  (4.50)  and  using  the  integration 

n 

result  for  cos  (a)da  ,  (4.38)  we  obtain: 


}  ' 


(1)0  a 1 


oh  (o)?  \ 


* 
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t  =  -  ~2  +  ~~  ^2  I  (P^  )(airn)p  cosP  nll£  sin  nll£ 
2^  k=o 
k  even 

2  k-  2m+l  ,  >-\  i  n/  \ 

n  •  (n-C)  +  0(e) 


m=l 


k-  2m+2 


-  +  C( n-5)  P1nn(a?n) 


1 

2 


,P+3 


INTC2^) 

P  y2  (P-1)  2k! 

k=o  2k)  (k!)2  22k 


(4.53) 


•  cosP  Pn]i5  sin  nllC  +  0(e) 


Also  t  (£,n)  =  -t_(n»?)  to  the  order  of  e  .  We  approximate  the  singular 

n  4 

curve  by  the  curve  J  =  0(e)  as  e  -*  0  ,  in  keeping  with  our  assumption 

oo  1c  1c 

that  the  sum  ^_|2et(5,n)  is  negligible.  So  the  singular  curves  are 
approximated  by 

t  =  0(e)  as  e  ->  0  and  t  =  0(e)  as  e  -*  0  . 

£  n 

For  q  =  1  ,  p  =  1  ,  we  have 

p  =  5  +  8[e(Ttn)2  a  sin  nll£]  1  .  (4.54) 

For  q  =  1  ,  p  =  2  , 

D  =  £  +  16[e(Hn)2  a2  sin  2nII£]  ^  .  (4.55) 

In  the  general  case  the  curves: 

_  1  +  P,qnn(ann)P  ■  0  ,  (4.56) 

2  2P+3 


where 


9ldx§xIgDn  i 
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INT(^) 

_  v  .P-1.  2k!  p-2k-l  _w 

M(5)  -  l  (  9v)  - 7~?t7  cos  nlK  sin  nn^ 

k=o  ZR  (k!)  2 

trace  out  the  approximation  to  the  set  of  points  in  the  (£,q)  plane  where 
J  =  0  .  If  we  choose  for  which  M(£  )  has  its  maximum  value  and  we 

choose  q^  such  that 

=  2P+2[epq(TIn)P+1aP  max  M(0  ]  1  ,  (4.57) 

we  have  an  approximation  to  the  minimum  time  for  which  J  =  0  . 

Increasing  or  decreasing  q  slightly  from  q^  ,  the  actual 
zero  of  J  and  t  should  occur,  as  (q-£^)  varies. 

Similarly  we  have  J  =  0  for 

q  -  £  =  2P+2[epq(Hn)i:)+^aP  max  M(q)]  ^  (4.58) 


when  t  =  0(e)  . 
q 

_  £- 

We  can  use  the  linear  result  for  t  ,  t  =  — —  to  obtain  an 
approximation  for  the  initial  time,  t^  for  which  breakdown  occurs. 

t,  =  2P+'*'[£pq(IIn)P+‘*'aP  max  M(q)]  ^  .  (4.59) 

b 

In  this  way,  we  arrive  at  a  prediction  of  initial  breakdown  time  for  the 
two  cases  studied  in  the  F.P.U.  report. 


For  p  =  1,  q=l,  n  =  1 ,  we  have 


, 


* 
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M(£)  =  sin  IT£ 

max  M(£  )  =  M(l/2)  =  1  (4.60) 


4 


This  result  is  in  agreement  with  the  papers  which  dealt  with 
this  problem.  For  p=2,  q  =  1  ,  n  =  1  ,  the  other  F.P.U.  problem, 

M(£)  =  sin  2nC 

max  M(£)=M(1/4)=1  (4.61) 


One  notes  further  from  (4.59),  that  initial  breakdown  time  t, 

b 

decreases  as  the  parameters  p,  q,  e,  n  and  a  are  decreased. 

Finally,  we  outline  the  calculation  of  Y(P)  at  a  point 

P(x,t)  .  We  calculate  £(P)  and  q(P)  from  (4.43)  and  find  the  point 

P  where  the  line  £  =  £(P)  intersects  x  =  0  or  t  =  0  .  When 
o 

£(P)  <  0  ,  q (P  )  is  determined  from  (4.42)  such  that  x (£ (P) , p (P  ) )  =  0  , 

and  we  have  Y(P  )  =  0  .  For  £ (P)  >  0  ,  P  is  on  t  =  0  and 

o  —  o 

YCPq)  =  a  sin  n£(P)  .  We  calculate  for  p  =  1,  q  =  1,  n  =  1. 

We  consider  the  first  term  of  these  series  expressions: 

all 

u(£,n)  =  ~  [c°s  Iln  +  cos  TI£ ] 

aTT 

v(£,n)  =  —  [cos  Tin  -  cos  TI£ ] 


(4.62) 
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t  (s»n)  =  +£ff“  (n-S)  sin  nn 

1  all^ 

x  (s,n)  =  -o  +£t t-  (n-C)  sin  nn 
n  z  ib 


(4.63) 


THE  (x  >  t )  PLANE 
Figure  XI 


We  assume  that  £(P)  <0  ,  and  we  calculate: 


Y(P)  = 


-  / 


n 


n(P) 

/  y,  dn 
n(P  ) 

o 

n  (P) 


[ux  +  vt  ]dn 
n  n 


n(P0} 

n(p)  n2 

J  [2  cos  n6  +  e  T7~  (6-?)  sin  116  ] d 6 


"<V 


(4.64) 
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For  small  t  and  small  e  ,  the  second  term  is  insignificant; 

for  t  =  0(— )  ,  n-£  can  be  considered  constant  over  the  interval 

e 

(n(P  ) ,n(P) )  »  so  for  Y(P  )  on  x  =  0  we  have 
o  o 


Y(P) 


—  sin  116 

n 


e  f-  (5-£(P))  cos  n6 

O 


n(P) 

+  0(e) 

n(po) 


(4.65) 


Similarly,  for  p  =  2  ,  we  have 


Y(P) 


—  sin  116 


(6-5(P))  cos  116 


n(P) 

+  0(e) 

n(Pc) 


(4.66) 
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CHAPTER  V 


THE  MOVING  THREADLINE 


We  now  attempt  to  apply  the  preceding  methods  to  the  problem 
of  determining  the  motion  of  the  'moving  threadline".  We  shall  investigate 
two  models  for  its  equation  of  motion;  for  details  of  the  derivation 
of  these  see  Zaiser  [18], 


2  2  _1 
7-  Y  +  aY  ^  *  (1  +  Y  )  Y  ,  Case  I  , 

4  XX  Xt  tt  X  XX 


2-1  -  —  -1 

7-  (1+Y2)  Y  +a(l+Y2)  2  Y  +Y^  =  (1+Y2)  ’  CaS6  11 

4  x  xx  x  xt  tt  x 


(5.1) 

(5.2) 


We  shall  impose  smooth  perturbations  in  the  boundary  conditions,  of  the 
form, 


Y(x,0)  =  s  sin  IIx 

Y  (x, 0)  =  0 
t 

Y(0,t)  =  Y(l, t)  =  0 


0  <  x  <  1 


t  >  0 


and  for  reasons  to  be  shown  later  we  take  0  <  a  <  2  . 


(5.3) 


(5.4) 


We  write  the  equations  as  a  first  order  system  by  setting 

Y  =  u  ,  Y  =  v  ,  and  thus  obtain: 
x  t 


/  0  ,1 

-12  — 
\  (l+u2)  -  f-  <1+U2)k  ,  -a(l+ii2)2 


(5.5) 


X 


" 


(  y  +  i)  *=  . . y  + 


Y*f  Y  * 

■ 
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where  k  =  0,  -1  in  Cases  I  and  II  respectively,  with  initial  conditions 


u  =  U(x,0)  =  e II  cos  ITx  , 
v  =  V(x,0)  =  0 


0  <  x  <  1  , 


(5.6) 


The  solution  for  the  "infinite"  threadline,  i.e.  ignoring  the 
equivalent  to  equation  (5.4),  can  be  obtained  in  a  manner  similar  to  the 
previous  problem.  But  the  end  conditions  for  the  finite  interval  cannot 
be  as  easily  satisfied  as  they  were  in  the  string  problem. 


We  look  at  the  corresponding  linear  problem  in  order  to 
establish  initial  conditions  for  the  infinite  string  which  will 
approximate  those  needed  to  provide  for  the  zero  end  conditions  in  the 
nonlinear  problems.  We  consider  the  equation 


Y..  +  aY  .  + 
tt  xt 


<f-  - 1) 


Y  »  0 

XX 


t  >  0  , 

0  <  x  <  1 


(5.7) 


with 

Y(0 , t)  =  g (x)  ,  0  <_  x  <_  1  ,  (5.8) 


where  g(0)  =  g(l)  =  0  , 


Yt(x,0)  =  0 


Y(0,t)  =  Y (1 , t)  =  0  ,  t  >  0 


(5.9) 


Setting 


u  = 


Y 

x 


V 

n 


and  v  =  Yfc  , 

+  <f  +  1)  uc 
+  <f  -  1)  -U 

2  P 


as  in  the  previous  problem  we  obtain: 

=  0  u(a,a)  =  f(a)  , 

(5.10) 


0 


v(a,a)  =  0 
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x5  -  <!  -  « 

h =  0 

x(a,a)  =  a  , 

(5.11) 

x  -  (f  +  1) 

rt- 

ll 

O 

t(a,a)  =  0  , 

n  2 

n 

where 

f  (a) 

will  be  defined  later,  and 

is  assumed  to  be  sufficiently 

smooth. 

We 

choose  Riemann 

invariants  r 

and  s  such  that: 

1 

1 

Sv  =  2 

r 

V 

2 

(5.12) 

3  =  J  (f  + 

u  2  2 

1)  r 

u 

=  _  I  _  u 

2l2  * 

Integrating  these  we  obtain: 

r  =  r  (£)  =  (1  -  ^)u  -  ^  , 

s  =  s(n)  =  i;  (1  +  ^  , 

so  that,  in  the  usual  way, 


(5.13) 


u(£  ,n)  =  r (C)  +  s  (n )  > 

v(£,n)  =  (1  -  -|)s(n)  -  (l  +  -|)r(5) 


(5.14) 


and  we  have,  from  the  initial  line 

r(e)  =  \  (1  -  f)f(0  , 

s(n)  =  ^  (l  +  f)f (n) 


(5.15) 


Again  from  (5.11)  we  have 

X  =  2  (1  +  2)n  +  2  (1  ~  2K 
t  =  \  (n  -  5)  ; 


(5.16) 


and  thus 


J  *jx  -nl  rnisinsJfcfl  saoo 


S 


<  (n)s  +  (;)*  =  Cn«5)u 

« 

.  <n)5(?-  +  I)  i  =  (n)e 
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n  =  x  +  (l  -  -|)t 
£  =  x-(l+-|)t 


(5.17) 


Now,  in  order  to  satisfy  the  end  conditions,  the  initial  conditions  are 
extended  to  (-00  ,  00)  in  the  following  fashion.  Let  h(x)  =  ,  anc} 
let 


h(x)  =  h(x)  ,  0  <_  x  <_  1  , 

=  hC^Y^)  ,  l<x  <_  1+y  (5.18) 

=  h(-yx)  ,  -  —  <  x  <  0  , 

Y  -  — 

OL  QL 

where  y  =  (1  -  ~rf)  /  (1  +  -j)  .  Then  the  initial  conditions  on  [0,1+y] 
are  extended  periodically  to  the  right  and  those  of  [-  y,  1]  , 
periodically  to  the  left. 


LINEAR  CHARACTERISTIC  LINES  FOR  MOVING  THREADLINE 


Figure  XII 


(x)rf  -  (x)rf 
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With  this  extension  we  see  that,  for  example,  at  the  point 
we  have  £  =  —  ,  p  =  1  ;  for  which 


r(-  i)  =  -  f)f(-  i)  -  -|(1  -  f)h(l) 

s(l)  =  |(1  +  |)£(1)  =  i(l  +  f)h(l)  , 


(5.20) 


and  thus 


v  =  |(1  -  f)(l  +  f)h(l)  -  |(1  +  j)(l  -  f)h(l)  , 
=  0 


(5.21) 


If  v  =  0  at  every  point  on  x  =  0  and  x  =  1  ,  then  we  have 

Y  =  /t  vdt  =  0  there, 
o 

Physically,  it  is  clear  that  a  <  2  .  It  is  also  evident 
from  the  characteristics  that  for  a  >_  2  the  problem  is  overdetermined . 

The  above  extension  does  not  provide  for  the  end  conditions, 
(5.9),  in  the  nonlinear  problems.  But  we  assume  that,  in  the  case  of 
small  perturbations  from  the  linear  problem,  it  approximates  the 
required  conditions. 


We  consider  the  nonlinear  problems  (5.6).  We  set  h(x)  = 
e II  cos  nx  ,  and  define 

k  _  JL 

,+  ct  / -|  .  2.2  2  2 

0=—  (1  +  u  )  +  (1  +  u  )  , 

k  _  1  (5.22) 

$  =  —  (1  +  u)  -  (1  +  u  )  , 

k  =  0,  -1  . 


x  tns 
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Then,  the  eigenvalues  of  the  matrix  in  (5.5)  are  , 

+  — 

and  left  eigenvectors  are  ($  ,1)  and  (f>  ,1)  respectively.  A  brief 

calculation,  as  in  previous  chapters,  provides  us  with  the  equations 
in  characteristic  coordinates, 


v^  +  $  u^  =  0 
v  +  0  u  =0 

n  n 


(5.23) 


-  $  t  =  0 
x  -  $+t  =  0 

n  n 


(5.24) 


The  initial  conditions  are  then  given  by 


u(a,a)  =  f(ar) 
v(a,a)  =  0  , 


(5.25) 


x(a,a)  =  a 
t(a,o)  =  0 


(5.26) 


We  choose  the  Riemann  invariants, 


r  =  r(£)  =  - 


v 

2 


u 


k 

2,2 


1 

2,  2- 


-  ±  !  [ja+rY  -  (l+yp  hdy 

o 


=  s(n)  = 


v  1  fU  .a..  2,2  ,  2  2,  , 

~2  +  ~2  J  [^(l+Y  )  +  (1+Y  )  ]dy 

o 


(5.27) 


so  that  we  have 


u 


r  +  s 


/  (1+y2)  2dy 


=  sinh  Xu 


(5.27) 
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s  -  r  =  v 


_k 

+  -^  /  (1+y2)2  dy 


,  a 

v  +-  U 


•  a  .  -1 

v  +  —  sinh  u 


k 

k 


0  , 

1  . 


(5. 28) 


Since  sinh  is  a  mono  tonic  increasing  function  of  u  , 

we  may  invert  the  equation  so  that: 


u  =  sinh  (r+s) 


/  Ci 

s  -  r  -  —  sinh  (r+s) 


v  =  < 


a 


v 


S  -  r  -  —  (r+s) 


k  =  0  , 

k  =  1  . 


(5.29) 


Again,  because  r  and  s  are  constant  along  E,  =  const,  and 
q  =  const.  respectively,  we  find  from  the  initial  conditions  that: 


and 


r(5)  =  -  f(?)  +  -j  sinh  1  f  (6) 


s(n)  =  ^  f  (q)  +  ^  sinh  f(n)  >  if  k  =  0 


(5.30) 


r(0  =  \  (1  ~  -|)  sinh  1  f(?)  , 


s(n)  =  (1  +  ^ )  sinh  1  f(n)  ,  if  k  -  -1 


(5.31) 


Thus,  the  solution  to  the  first  part  of  the  problem  can  be 


expressed  in  the  form: 


:rt£rto  bucj:..  . 
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u(?,n)  =  sinh  [-|(f  (n) -f  (5))  + -|  (sinh  1f(£)  +  sinh  1f  (r>)  ] 

(5.32) 

v(?,n)  =  -|(f  (S)+f  (n))  -  (sinh  1f(n)  -  Sinh  1f(?))--|  u(?,q) 

for  k  =  0 

and 

u(?,n)  =  sinh  [-|(1  -  ^-) sinh  ^f(?)  +  ■j(l  +  y)sinh  1f(n)] 

(5.33) 

v(? ,n)  =  4(1  -  -f)  (1  +4^  [sinh  lf  (n)-sinh  1f (?)  ]  for  k  =  -1  . 


with 


We  now  turn  to  consideration  to  the  second  system  of  equations 


-  $  t^  =  0 


x  -  $+t  =  0 

n  n 


(5.34) 


x(a,a)  =  a  t (a, a)  =  0 


(5.35) 


2  2 

Firstly  we  note  that  (1+u  )  =  cosh  (r+s)  so  that 


$  =  [cosh  (r+s)]^  +  [cosh  (r+s)]  ^ 


(5.36) 


and 


r(?)+s(n)  =  -|  [  f  (n ) -f  ( S )  ]  +  -|[sinh  1f(?)+sinh  1f(n)]  , 


r(?)+s(n)  =  -|(1  -  y)  sinh  1f  (?)  +  -|(1  +  -|)sinh  1f(n)  , 


if  k  =  0  , 
(5.37) 
if  k  =  1 


+ 


In  order  to  express  $  ,  $  in  series,  we  set  f (x)  =  em(x)  , 

and  note  that  if  we  maintain  e  <  ,  so  that  | r+s |  <  ^  ,  then 


0  «  J  * 


8+i )  riaoo  **  (  u+I) 


51  ,  [Crr)li“rinxB+(3)l''rinjaJ7  +  [<3)i-(ri)l]  ■  -  (fc)e+(3)-i 
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+ 

cosh  (r+s)  is  an  analytic  function  of  (r+s)  and  the  series  for  $ 
is  convergent. 


r(£)  +  s(n) 

[cosh  (r+s) ]  ^  = 
Thus  we  have 


=  |(1  +  |)m(n)  +  -|(1  -  f)m(C)  +  0(c3)  , 

2  9  / 
!-—[(!+  f)m(n)  +  (l  -  f)m(0]  +  0(e). 


(5.38) 


(5.39) 


+  2  , 

$  =  f±  (1  -  f-  c(£,n)  +  0(e  ))  if  k  =  0 


$  =  (|  +  1)  (1  -  C(C,n)  +  0(e'  ))  ,  if  k  =  -1 


(5.40) 


where  we  define 

C(?,n)  =  [(1  +  f)rn(n)  +  U  -  f)m(5)]2  .  (5.41) 

+ 

Proceeding  as  before  with  P  =  $  and  Q  =  $  ,  we  integrate  (5.34) 

to  find 


x-n  =  (^  -  l)t(c,n)  +  -|—  (1  -  f )  k 

2  V 

X-C  =  (*2  +  l)t(£r,a)  -  g-C-j  +  1) 

and,  on  subtracting,  we  have, 

_r  2  n  , 

t(5,n)  =  2  "16  ^  [(1  "  c(°>n)to(a,n)-(f  +  1)  KC(S ,o) t0 (?,a) ]do 

+  . . .  ,  k  =  0,-1  .  (5.43) 


5 

/  C(a,n)ta(a,n)da  +  ... 
n 

(5.42) 

n  :  1 

/  C(£,cr)  t^(£,cr)da  +  . . .  , 


We  look  for  solutions  of  the  form: 
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t(5>n)  =  l  e21  (£,n)  • 


L 

1=0 


Taking  the  first  two  terms  of  this  series  we  obtain 


2  n 

t(;,n)  =^  +  -§2  /  I(f  +  1)  C(c,0)+(1  -  -f)  KC(cr,n)]dcr  ,  (5.44) 

x(s,n)  -  ^  +  |  t(5,n)  +  -§2  Ad  ”  |)"kC(a,n)-(f  +  l)~kC(S  a)]do,  (5. 


and  so  we  may  write 


x(? , n)  =  -|(1  +  |)n  +  ■§(!  -  f)C 


2  n 


+ 


32 


/  tcf 


l)(f  +  D"kC(S,a)  +  (|  +  1)  Cl  -  |)  kC(a,ri)]d0 


k  =  0,1 


(5.46) 


3  /  x  t ) 

If  the  Jacobean  J  =  — j — r-  ^  0  ,  we  may  solve  for  £(x,t)  , 

o  fs  5  h ) 

n(x, t)  .  Let  us  seek  expressions  of  the  form 


£(x,t)  =  £  (x,t)  +  e  £^(x,t)  +  ... 


n(x,t)  =  nD(x,t)  +  e  n1(x,t)  +  .. 


•  • 


(5.47) 


We  then  find  £  =  x  -  (1  +  ybt  ,  n  =  x  +  (1  -  Ar)  t  .  Combining  (5.44) 

o  2  o  z 

and  (5.45)  we  have 


5  =  x-(l  +  §)t-e2[x(1)U,n)-(l  +  |)t(1) (5,n) ]+. • • 
n  =  x+(l  -  |)t+c2[x(1)(5,n)+(l  -  § )t (1) (5 ,n) ]+. • • 


(5.48) 


45) 


From  this,  we  have 


sixto  io  aura  <3 
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2  x+(l  -  -|)t 

S  =  S(x,t)  =  x-(l  +  f)t  +  jt  /  (-  +  l)“kC(?,a)da+. ..  , 

x-(l  +  |)t  2 

(5.49) 

2  x+(l  -  -|)t 

n  =  n(x,  t)  =  x+(l  -  —) t  +  ~  /  (1  -  —)  kc (a, n)da+. . . 

x-(l  +  -|)t 

Finally,  we  have  to  solve,  say  by  successive  approximations, 

2 

£  =  Ux,t)  =  x-(l  +  |)t  +  F(S)  +  0(e2) 

2  (5.50) 

n  =  n(x,t)  =  x+(l  -  -|)t  +  E(n)  +  o(e2) 

where 

"k  2  ?  2  v2 

E(n)  =  (1  '  f)  Id  +f)"»(n)  +d  -f)  y"\  , 

-k  2  ,  2  2  (5-51) 

F()-(l+|)  [  (1  -  f)  m(C)  +(1  +  f)  y"\  • 

It  is  to  be  expected  that,  in  the  neighbourhood  of  t.  there 

b 

will  be  difficulty  in  solving  expressions  such  as  (5.50)  because  p  p 

X  t 

are  unbounded  there. 


Thus  we  have  £(x,t)  ,  p(x,t)  ,  and  combining  this  with  the 
results  u(£,p),  v(£5p)  ,  (5.32),  (5.33)  we  have  u(x,t)  ,  v(x,t)  , 

which  one  may  integrate  to  obtain  Y(x,t)  . 

These  results  are  based  on  perturbations  on  the  linear  problem. 
If  t  is  small  and  e  is  small,  clearly  the  approximation  is  good. 

However,  if  t  is  of  the  same  order  of  magnitude  as  — ,  we  must 

c 

examine  to  what  extent  boundary  conditions  are  satisfied. 
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In  the  linear  case,  the  line  x(£,n)  =  0  has  the  equation 


(1  +  |)n  +  (1  -  |)5  =  0  . 


(5.52) 


But  (5*46)  shows  that  these  values  are  altered  by  a  term  which 
for  large  values  of  (q-£)  could  have  considerable  effect.  If  we 


consider  a  point  P(5^»n^)  lying  at  the  intersection  of  (5.52)  wi 


th 


xC 

q  -  £  =  — -  ,  where  k  is  chosen  depending  on  the  value  of  t,  which 
Z  b 

£ 

we  consider  later,  we  find  the  deviation  of  x(£^,ri^)  from  x  =  0  . 
With  the  condition  (5.52),  and  m(x)  taken  from  equation  (5.10),  and 
the  fact  that  f  (x)  =  etn(x)  ,  we  integrate  (5.46)  to  obtain: 


,a 


xUvr,1)  =  -  1)[(1  -  f)2m(e1)2  +  f-(l  +  -f)2] 


+  (f  +  1)  [  (1  +  f)2m(ni)2  +  f-(l 


f)2]}  +  0(E2) 


k  =  0 


(5.53) 


We  note  that  f(n^)  =  f(£^)  along  (5.52)  and  take 


,  =  _ £_ 

23 

e  11  a  +  f) 


(5.54) 


so  that 


I  ,r  \  |  (2-a) (a  +10a)  .  nn 

|x(c1,n1)|  l - o -  <  -08 

1  1  2h(2+a) 


0  <  a  <  2 


(5.55) 


Similarly,  for  k  =  -1  ,  we  set  (n^“£^)  =  "2  3 


16 


e  (i  +  ^)3 


and  we 


- — 
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calculate  x^i,r,l^  <  *05  for  0  <  a  <  2  . 


Thus  it  seems  we  have  a  reasonably  close  approximation  to  the 


given  boundary  line.  The  analysis  for  x  =  1  is  similar.  We  note 
further  that  along  (5.52)  we  have  f(£)  =  f(p)  ,  so  that  from  the 
expressions  for  v(£,p)  in  (5.33)  and  (5.34),  we  have  v  =  0  along 
that  line. 


Thus  if  we  accept  the  first  approximation  to  x  =  0  ,  we  see 


that  our  boundary  conditions  are  satisfied  exactly.  If,  on  the  other 
hand,  we  wish  to  proceed  to  the  second  approximation  to  x  =  0  ,  this 
will  not  be  so.  In  this  case,  by  the  continuity  of  v  ,  it  is  clear 
that  v  will  not  differ  appreciably  from  zero.  However  the  value  of 
Y  may.  This  point  would  appear  to  require  further  study. 


Now,  we  have  breakdown  if 


J  =  [P-Q]tr t  =  0 

£  n 


(5.56) 


From  (5.44)  we  have 


(5.57) 


where 


C  (a,n)  =  2(1  +  yp)  [  (1  +  yp)m(n)  +  (1  +  ^)m(a)  Jm’  (p)  .  (5.58) 

p  l  L  Z 


Thus 


t  =  -  \  +  fjd  -  f)  kd  +  f)2n>(n)m'(n)(n-0  +  o(e2)  , 


(5.59) 


ii  .tv  •>[££•  id  >  T  .  s  «v 
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and  we  have  =  0  +  0(e  )  when 


n-5  -  [d  -  f)'kU  +  -|)2m(n)m'(n)]'1 


(5.60) 


The  factor  m’  (n)  takes  on  three  forms  which  are  given  by  (5.18).  In 
order  to  obtain  the  maximum  value  for  m(n)m!  (ri)  ,  we  consider  the  form 


m(n)  -  n  cos  (n  -i^)~rl) 


(5.61) 


j ( 1+Y ) ~ Y  <  n  <  j (1+y)  j  =  1,2,... 


so  that 


m(n)m’  (n)  =  sin  ( 211  j  n) 


£ 

-  2y 


(5.62) 


o  n 

We  employ  the  approximation  t  =  t  =  — y-  to  determine  the 

time  of  breakdown  t^  ,  and  we  take  the  minimum  value  of  satisfying 

(5. 60) ;  thus 


ct  “ 


8(1  -  f)1+k 

2  3,,  ,  a, 3 

e  n  (1+y) 


k  =  0,-1 


(5.63) 


A  consideration  of  the  term  t^^  in  (5.44)  shows  that 


e2  t(1)l  <  2 


(5.64) 


for  time  t,  ,  so  that  where  e  is  small,  t,  =  0(— r-)  ,  and  (5.63) 

b  b  2 

e 

is  significantly  representative  when  breakdown  occurs. 


Finally,  we  indicate  a  procedure  for  determining  the  solution 


Y(x,t)  to  the  problem.  Consider  a  point  P(x,t)  where  0  <  t  <  t^  and 
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n(PQ)  =  -  yC(P) 


(5.68) 


Figure  XIII 


Now,  the  range  of  integration  of  (5.65),  from  n(PQ)  to  n(P) 
is  not  large,  so  it  suits  us  to  consider  only  the  first  terms  of  each 
series  expression.  We  perform  the  sample  solution  for  the  first 
threadline  problem,  k  =  0  ; 

u  =  f[(1  +  f)m(p)  +  (1  -  f)m(0] 

2 

v  =  -|(1  -  |-)  [m(n)-m(n)  ] 

1  2 

x  =  4(1  +  4)  +  TT  (n-S)  (4  +  1)  m(n)mT  (p) 

p  2  2  16  2 

1  2  o 

t  =  4  +  TT-Cn-OOf  +  !)  m(n)m'  (p) 

p  2  16  2 


(5.69) 
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We  set  Y(Pq)  =  0  because  of  the  desired  boundary  conditions 


and  thus 


Y  = 


n  (p) 

/  -^(1  +  ^)m(6)d<5 

n(P  ) 

o 


(5.70) 


n(p)  3  , 

+  /  “(6-0  d  +^)  [(l+Y)m/(6)m,(6)+2Yin(C(P))m(6)in'(6)]d6 

n(P  ) 
o 

Clearly,  for  t  <  ,  when  e  is  small,  the  second  term 

is  insignificant. 

The  calculations  are  straight  forward  except  for  the  fact 
that  m(x)  is  defined  differently  on  different  intervals;  that  is,  from 
(5.18) 


m(x)  =  II  cos  n  (x-x . )  ,  x  .e  I.  ,  i  =  3,4,7,8,11,12,... 

i  l 


=  n  cos  n 


X  ,-x 
1 

Y 


x  £  I.  ,  i  =  2,6,10,. . .  (5.71) 


=  II  cos  iry(x-x.)  ,  x  £  I.  ,  i  =  1,5,9,... 
1  1 


Let  us  suppose  that,  as  in  Figure  XIII,  n(P)  £  IQ  ,  n(P  )  e  I. 

o  o  4 

so  that  we  must  consider  the  integral  for  each  interval  1^,  1^,  Ig  . 


From  (5.67)  and  (5.68), 


n(P)  =  x  +  (1  - 


-)t 

V 


(5.72) 


n(PQ)  =  -y(x 


a. 


<1  +  f)t 


and  from  (5. 70) , 


ab(S)m(£  +  «!■ 


76  - 


x. 


x 


8 


y  =  f  n  cos  n(x-x.)dx  +  J  n  cos  n 


n(P  ) 

o 


X , 


X  ,-x 
6 

Y  ' 


n(P) 

+  /  n  cos  n(x-x8)  , 

x8 


where 


Xi  =  ’  i  =  4,8,  . .  .  , 

(  r\ 

xi  =  1  +  — —  (1+y)  ,  i  =  2, 6, 10, .  .  . 

As 

sin  II  x±  =  0  ,  i  =  1,2,3,... 

there  remains  only 

Y  =  sin  II(x+(l  -  ^)t-x^) 

-  sin  II  (~y  (x- (1  +  -|)t)-x8) 


(5.73) 


(5.74) 


(5.75) 


In  general  then,  computations  involve  a  consideration  of  the 
particular  intervals  over  which  calculations  are  made. 


CHAPTER  VI 


NUMERICAL  RESULTS 


The  object  of  the  present  chapter  is  to  calculate  and  plot  various 
members  of  the  families  of  characteristic  lines.  The  initial  intersection 
of  two  members  of  the  same  family  of  lines  indicates  that  t  =  0  or 

t  =  0  at  the  point,  thus  we  have  J  =  0  and  the  occurrence  of  break- 

in 

down.  For  the  fixed  string  problem,  we  have  plotted  a  few  examples  to 
indicate  the  path  of  the  lines  and  the  shape  of  the  singular  curves  where 
J  =  0  ,  and  we  have  compared  the  breakdown  times  obtained  in  this  way, 
to  the  analytical  results  of  Chapter  IV. 


Of  the  methods  available,  the  method  of  Masseau  [8]  was  selected 
for  its  simplicity,  and  for  its  accuracy  in  view  of  the  slow  variation 
in  the  slope  of  the  characteristic  lines.  The  method  consists  in  a 
calculation  of  points  where  characteristic  lines  intersect,  and  thus 
provide  an  approximation  to  the  path  of  characteristic  lines  by  joining 
the  points. 

Consider  a  system  of  two  first  order  equations  in  two  dependent 
and  two  independent  variables  (2.1).  Having  written  it  in  terms  of 
characteristic  parameters  (2.5),  (2.6)  we  replace  the  derivatives  with 
expressions  such  as: 


us  = 


u(P1)  -  u(P2) 


P  “  P 
1  2 


9 


(6.1) 
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and  thus  obtain 


u(P3)  -  u(P2) 

Tp  -  p  I 

1  3 


+  b 


v(P3)  -  v(P2) 

lP3  -  P2l 


=  0 


u(P3)-  u(P1) 

lP,  "  Pi I 


+  d 


v(P3)  -  v(p  ) 

fl  -  p  I 

1  3  l1 


=  0 


X(P3}  "  X(P2}  t(P3}  “  tfP2} 

+  e  - m - rr— i -  =  0 


P  -  P 
3  2 


Pn  -  P^ 
3  2 


x(P  )  -  x(P  )  t(P  )  -  t(P  ) 

+  f  - ri - z  r~  =  0  , 


P  -  P 
3  1 


P  -  P 
3  1 


(6.2) 


(6.3) 


where  a,  b,  c,  d,  e,  and  f  are  functions  of  u,  v,  x,  and  t  , 
pi»  P2>  and  P^  are  ajacent  grid  points  as  in  Figure  XIV  and  |  p_^  -  p 
is  the  distance  between  the  two  points. 


Figure  XIV 
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Consider  two  points  P^,  (such  as  points  of  the  initial 

line),  where  u,  v,  x,  t  are  known.  Calculating  a,  b  and  e  at  p2  , 
and  c,  d  and  f  at  P^  ,  we  solve  (6.2)  and  (6.3)  for  u(P3)  >  v(p^), 
x(P^),  and  t(P^)  •  This  serves  as  a  first  approximation  to  the  values 
for  the  point  of  intersection  of  the  lines  through  P1  and  P2 


Then  we  take 


e  =  \  [e(P3)  +  e(P2)] 

f  =  \  [f (P3)  +  f(P1)] 


(6.4) 


and  by  (6.3)  calculate  new  values  for  x(P3),  t(P3)  .  Using  similar 

averaging  for  a,  b,  c,  and  d  which  are  calculated  with  the  most  recent 
values  of  u,  v,  x,  and  t  ,  we  solve  for  u(P3)  »  v(P^)  from  (6.2). 
Alternating  these  last  two  steps,  the  solution  of  x(P3),  t(P3),  and 
of  u(Pg)  >  v(P3)  »  should  give  a  convergent  series  for  u,  v,  x,  t 
at  P3  . 


Proceeding  from  a  set  of  points  on  the  initial  line,  a  grid 
of  points  is  built  up  approximating  the  points  of  intersection  of  the 
characteristic  lines.  (Figure  XIV). 


In  the  problems  we  deal  with,  we  can  determine  u  =  u(r+s)  = 
u(r(£)  +  s(n) )  from  the  Riemann  invariants  corresponding  to  the 
characteristic  lines  which  intersect  at  the  point,  and  we  can  determine 
r(£)  ,  s(n)  from  initial  conditions.  So  it  is  unnecessary  to  employ 
the  above  iterative  process  in  order  to  determine  u  ,  and  thus  the 


(eqa:Je  ow3  Jsfil  Saadi  gniJanm)  A 
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slopes  -e  ,  -f  of  the  characteristic  lines,  which  are  functions 
of  u  only,  are  known.  If  we  take  the  average  slope  (6.4),  of  the 
two  end  points  of  the  line  segment,  we  can  calculate  x(P^),  t(P^) 

immediately. 


The  inversion  of 

r  +  s  =  /  (1+  euP)^^2  ,  (6.5) 

o 

to  provide  u  =  u(r+s)  cannot  be  carried  out  exactly,  except  for  p  =  1.  But 
for  the  finite  number  of  discrete  values  r  and  s  which  enter  into  the 
calculations,  we  can  use  Einstein's  method  to  approximate  a  value  of 
u  for  which 


f(u)  =  r  +  s 


(6.6) 


where 


p+1  2  2p+l 

f  (u)  =  u  +  ^  e  +  £(£_n  + 

UUJ  u  2  p+1  2  4  }  2!  •** 


(6.7) 


The  Einstein  approximation  is  simple.  As 


(r+s)^  =  (1  +  eu^)^2 


u  < 


(6.8) 


the  function  r+s  is  monotonic  in  u  .  Incrementing  u  by 


„  r+s  -  f (u) 

Au  =  - jx  , 

(1  +  eu  ) 


(6.9) 


gives  a  convergent  sequence  of  values  of  u  .  The  degree  of  accuracy 


(u3  +  £) 
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of  the  numerical  results  is  increased  by  taking  a  denser  set  of 
characteristic  lines,  and  thus  shortening  the  straight  line  segements. 
As  the  slopes  of  the  lines  vary  slowly,  the  nonlinearity  being  a  factor 
of  e  ,  the  straight  line  segments  with  mean  slope  are  a  good 
approximation  to  the  curve.  An  indication  of  the  accuracy  of  the 
numerical  results  was  found  in  a  comparison  of  the  results  obtained 
when  the  number  of  initial  points  in  0  _<  x  <_  2  were  given  the  two 
values  64  and  80  .  The  resulting  times  of  breakdown  differed  by 
at  most  .03  ,  and  corresponding  x  coordiantes  differed  by  at  most 
.003  .  Thus  the  numerical  results  would  appear  to  be  satisfactory. 

The  following  is  a  comparison  of  the  values  of  £,  p,  t  at 
initial  breakdown,  as  arrived  at  by  the  analysis  of  Chapter  IV,  and 
the  numerical  work  on  the  fixed  string  problem. 

For  q  =  1  and  n  =  1 


II 

h-1 

Analytical 

Numerical 

£ 

a 

£ 

n 

t 

£ 

n 

t 

.3 

1 

.5 

3.2 

1.45 

.54 

2.  58 

1.09 

.1 

1 

.5 

8.62 

4.06 

.50 

8.50 

3.97 

.01 

2 

.5 

41.10 

20.50 

.50 

42.50 

20.  90 

p  =  2 

.3 

1 

.  25 

1.97 

.86 

.15 

2.62 

.91 

.1 

1 

.25 

5.41 

2.38 

.18 

5.26 

2.67 

.01 

2 

.  25 

13.15 

6.45 

.25 

16.10 

7.64 

TABLE 


I  =  n  bne  I  *  p 


• 
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As  is  seen  from  the  tabulated  results,  for  fairly  large  e  , 

when  for  p  =  1  ,  the  results  for  numerical  and  analytical  calculation 

of  breakdown  are  close,  except  where  e  is  very  close  to  - - - 

(anII)P 

For  p  =  2  the  linear  approximation  to  t^  is  representative,  but 
less  accurate.  The  two  graphs  (Figures  XV  and  XVI)  indicate  the  curves 
that  the  characteristic  lines  follow.  The  intersection  of  the  lines 
indicates  breakdown. 
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COMPUTER  PLOT  OF  CHARACTERISTIC  LINES  FOR  THE  FIXED  STRING 
e  =  .  3 ,  a  ®  1 ,  p  =  1 ,  q  “  1 ,  n  *  1 


Figure  XV 
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COMPUTER  PLOT  OF  CHARACTERISTIC  LINES  FOR  FIXED  STRING 
e  *=  o  3 ,  a  =  1,  p  =  2 ,  q  =  1,  n=l 


Figure  XVI 


CHAPTER  VII 


CONCLUSIONS 


It  appears  that  we  have  indicated  the  solutions  to  the  problems 
we  have  dealt  with,  but  there  remains  a  number  of  unanswered  questions. 

We  have  found  a  solution  to  (1.120  for  the  fixed  string,  and 

shown  that  after  time  t,  ,  the  second  derivatives  of  Y  ,  are  unbounded. 

b 

Thus  the  solution  cannot  be  extended  smoothly  beyond  such  points.  As 
a  result,  the  recurrence  phenomenon  observed  in  the  F.P.U.  report  [7] 
cannot  be  studied  by  means  of  the  model  we  have  studied.  Further  work 
in  this  has  been  carried  out  in  [16]  by  employing  a  different  model. 

We  have  examined  two  models  for  the  moving  threadline,  using 
the  initial  conditions  that  would  give  the  required  end  conditions  for 
the  unperturbed  problem  with  straight  line  characteristics.  We  have 
estimated  a  time  t^  at  which  breakdown  occurs  in  these  problems  as 
well.  But  as  we  have  shown,  the  lines  along  which  we  have  set  v  =  0 
deviate  somewhat  from  x  =  0  and  x  =  1  .  This  is  indicative  of  the 
fact  that  we  have  not  examined  precisely  the  boundary  value  problem  we 
originally  stated. 

An  improvement  on  our  approach  would  be,  as  we  stated,  to 
vary  the  initial  conditions  slightly  so  that  we  have  v  =  0  along 
precisely  x  =  0  and  x  =  1  .  This  would  involve  beginning  with 
the  interval  [0,1]  ,  and  determining  the  pairs  of  characteristic  lines 


II V  flariAHJ 
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which  intersect  at  the  boundaries  to  a  greater  degree  of  accuracy  than 
we  have  already.  Then  the  initial  condition  is  assigned  so  that  the 
Riemann  invariants  combine  as  required  to  give  v  =  0  .  Difficulties 
in  actually  carrying  out  these  calculations  have  prevented  us  from 
employing  this  procedure  and  the  matter  is  open  to  further  study. 

The  numerical  calculations  in  the  F.P.U.  report  [7]  and 
Zaiser's  Thesis  [18]  failed  to  indicate  the  singularities  we  predict 
in  our  analytical  and  numerical  results.  This  is  attributed  to  the 
nature  of  their  approaches  which  tend  to  smooth  out  discontinuities. 
Rather,  it  is  an  approach  such  as  ours,  which  examines  the  behaviour 
of  characteristic  lines,  which  is  successful. 


!c 

n3 qo  el  ie3±  rn  rfo  bna  arut  >cn  '  b  d  clvolci 

, 
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APPENDIX  I 


EQUIVALENCE  OF  THE  TWO  FORMS  OF  PERTURBATION 


Consider  the  two  methods  of  perturbation  we  have  used.  Firstly 
we  consider  problems  of  the  form 


V  =  F(eV£U 

u(x,0)  =  f(x) 
ut(x,0)  =  0 


t 


5 


)u 


XX 


(A. 1.1) 


and  secondly  of  the  form 


At  =  F(WVXX 


v(x,0)  =  ef(x) 


vt(x,0)  =  0 


(A.1.2) 


We  assume  that  for  all  values  of  x,  t,  e  considered  that 
solutions  do  exist,  and  are  unique.  Then  these  two  problems  are  equivalent 
in  the  sense  that  given  a  solution  of  one  of  them  then  the  solution  of 
the  other  may  be  found  by  the  relation 

v  =  e  u  ,  (A. 1.3) 


unless  of  course  e  =  0  .  In  this  case  it  is  obvious  that  v  =  0  and 

if  lim  u(e,x,t)  =  u(0,x,t)  then 
e  o 

=  lim  v(e,x,t)  =  lim  e  ti(e,x,t) 

e  -*•  o  e  o 


o 


.  .f  4  9i).t  lo 
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It  is  not  clear  however  that 


lim  u(e,x,t)  =  lim  —  v(e,x,t) 

£ 

£  ->  O  £  ->  O 


and  in  this  sense  the  equivalence  may  fail. 
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APPENDIX  II 


CONSIDERATION  OF  THE  SUCCESSIVE  APPROXIMATION  RESULTING  FROM  (2.37) 


Consider  equation  (2.37)  and  assume  that  we  have 


k-1  _(k) 


P=  l  £K"XP 

k=l 


k-1  „(k) 


Q=  l  e*-XQ 

k-1 


(A. 2.1) 


If  we  substitute  these  series  in  the  integral  equation  (2.38)  we  find 
that 


Tj  oo 

r-r  /  I  ek  [P(k)(a,n)t „(a,n)+Q(k)(C,o)trr(;,a)Jda  .  (A.2.2) 

p°-Q°  p°-Q°  k-i  a  0 

The  expression  for  t  in  [16]  is  a  power  series  in  e  .  If 
00  n  ( n ) 

we  set  t  =  e  tv  in  (A.2.2),  and  formally  rearrange  the  terms, 
we  obtain 


t  = 


n-g 

o  o 
P  -Q 


o  ^o 

p  -Q 


/ 


n 


I  s11 1 

n=l  k=l 


(k)  (n_k)  (k)  (n_k) 

[P  ;(a,n)ta  (a,n)+Q^  ;(S,a)ta 


(5»cr)  ]da 


(A.  2. 3) 

Provided  now  that  the  series  is  uniformly  convergent  for  values  of  £,p 
considered  we  can  interchange  the  order  of  summation  and  integration  to 
yield  (2.32).  The  convergence  of  series  concerned  is  considered  in 


Appendix  III. 


■ 
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APPENDIX  III 


CONVERGENCE  OF  THE  SERIES  FOR  t,  x,  t^,  t 

V  n 


We  will  show  that  the  series  for  t,  in  (2.31),  converges 
under  the  conditions  set  out  below.  Our  aim  is  not  to  find  the  most 
general  conditions.  For  simplicity  in  the  following  we  assume  _>  0  . 

First  we  note  that 


Jn  P(k)a,a)t^n)(£,cr)da  =  p(k)(£,n)t(n)(£,n)-  f  P^k)  (5,a)t(n)  (S,a)da 


(A.  3.1) 


We  begin  with  the  expression  given  by 


p  k 

[P 

K  n=l  a 


t<k)  =  [P°-Q°]_1  <!  /"  !  [P(n)(E.a)t(tn>(5.a)-Hl(n)(o,n)t(^n)(a.n)Ma 


-  I  [P(n)(«,n)-Q<n)(C,ri)]t(k  n)(5.n) 


(A.  3. 2) 


n=l 


N  ^  (^) 

We  define  S„  =  ,  Z_  e  |t  j  ,  and  note  that  for  given  quantities 


N  k=o 


a^k^  ,  that 


V  _ (k) | fc(n~k) 


I  L 

n=l  k^l 


e  a  t 


l  l  ena(k)|t(n’W 
k=l  n“k 

f  ck.(k) 

k-1 


e  a  S 


N-k 


(A. 3. 3) 


Now  consider  the  sum  S.T  .  We  have 

N 
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S._  <  .  |  ' 

N  -  p°_Q°  | p°_Q° 


N  k 


I  I  Ek|p(n,(5,n)-Q(n,(5,n)||t(k-n>(5,n)| 


k=l  n=l 


l  ek  /n  I  P^n)(5,a)t(k  n)(5,a)+Q(n)(a,n)t(k  n)(a, n)  do 
k=l  S  n=l  a  a 


n-C 


N  N  .  ,  , 

n  V  k-n.  (n) 


P°-Q°|  |P0-Q° 


l  en  l  eM|pw(5,n)-Qw(C,n)l|t 


(n)  tr  «<vl  [  *.  (k-n) 


n=l  k*n 


N  N 

l  ^ l 

n=  1  k=n 


I  cnj?  ek-nf  |P^n)(5,a)t<k-n)a)o)Vn)(o,n)t(k-n)(a,n)|da 

5  a 


<  s  + - - -  l  en|p(n) (5,n)-Q(n) (S,n)  | s 


°  JP°-Q°!  n=l 


N-n 


N  n 


-i—  I  /  en|P<“>«,a)|sH.n+|Q(»)(0,n)|s  i..  (A.3.4) 

P  -Q  I  n=l  £  v  u 


is  monatonically  increasing  of  N  so  that 


S  <  S  -  ,  n  >  1 

N-n  —  N-l  — 


We  further  stipulate  that  we  are  considering  those  values  of  e  for  which 
the  series  for  P  and  Q  are  absolutely  convergent,  and  assume  a  range 
of  values  of  £,r|  such  that 


V  k  ip(k)| 

2  e  |P  I  i  A  , 

k-1 


2  e  |Q  I  1  B  , 

k-1 


oo 

I  ek|P.(k)|<A'  , 

k-1  ^ 


(A.  3. 5) 


I  sk  lQr(k)l  <  B'  . 


k-1 


‘5  1  ~ 
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Thus 


S.T  <  S  + - - -  (A+B)  S  _  +  j 

N-  °  |po_Q0|  j  N-l  J 


(A'SN-l+B'SN-l)d° 


(A.  3.  6) 


It  remains  to  show  that  there  exists  a  K  such  that 


SN£K  e 


n-5 


We  do  this  inductively.  First  we  note  that: 


1  <  en~g 

S°  "  |P°-Q°i  "  |P°-Q° 


(A. 3. 7) 


Then  we  assume  that  K  >  |P°-Q°|  ^  >  0  ,  and  further  that 


S  <  K  e 

n 


n-5 


n  <  N-l 


From  (A. 3. 6)  we  have 


_  i _o  _o | — 1 

SN  <  |P  -Q  I 


n 


{(n-5)  +  (A+B)K  e11  S 


+  j  A1  K  e°  B'  K  e"  ^da} 


n-5. 


(A. 3. 8) 


and  clearly 


SXT  <  I  P°-Q°  i  1  { 1+K  (A+B+A ’  +B  ’  )  en  . 

N  —  1  1 


(A. 3. 9) 


We  require  a  K  >  0  such  that 


1  +  K (A+B+A '  +B  ’  )  <  K  |P°-Q° 


K  _>  [  j P°-Q°  |  -  (A+B+A'+B')]  1 


(A. 3. 10) 
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Thus,  the  required  K  exists  under  the  condition  that 


o 


(A+b+A’+B')  <  | P  — Q 


(A.  3. 


Clearly,  we  can  choose  e  small  enough  so  that  this  holds. 

Thus,  provided  n~£  is  bounded,  the  series  is  bounded.  As 
is  monotonic  and  bounded,  it  must  converge  to  a  finite  limit,  for 
each  value  of  (£,n)  •  We  can  show  in  a  similar  way  that  the  series 


for  x,t^  and  t  converge  under  the  same  conditions. 

£  n 


ID 
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APPENDIX  IV 


ORDER  OF  OMITTED  TERMS  IN  SETTING  t  =  At*'1'* 

£  £  £  r 


Our  assumption  in  predicting  the  time  of  breakdown  is  that 

00 

£  £ k  t^k)  <  e  K  (A. 4.1) 

k=2  ^ 

for  some  constant  K  .  If  a  specific  representation  of  t(£,p)  is  known, 
as  in  Zabusky  [15],  it  is  possible  to  check  the  assumption  (3.7).  We 
shall  attempt  to  verify  our  procedure  here  with  respect  to  the  problem 


y 

y 


tt  -  a  +  Eyx)y: 

(x,0)  =  a  sin  II  x 


y  (0 ,  t)  =  y  (1 ,  t)  =  0 


(A. 4. 2) 


If  we  reduce  this  problem  to  a  set  of  equations  as  in  Chapter  II,  it 
follows  that,  if  we  set 

1 

f(x)  =  (1  +  call  cos  Ex)2  ,  (A.  4.3) 

then  we  find  the  equation 

f  ,  f2(n)  f*(n)  t  ,  f2(£)  f 1  (£)  ,  =n 

nC  f3(0  +  f 3  (ft)  ?  f3(£)  +  f3(n)  n 

t(£,p  =  0 

te(£,£)  =  -  tn(£,0  -  -  \  f(£) 


(A. 4. 4) 


5  * 

' 
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Again  if  we  write  t  =  T  +  0  with 


T(5,n) 


-?/ 


n 


dy 


f(y) 


(A. 4. 5) 


then 

T  =0 

5n 

T(5,5)  =  0 

Tr(5,5)  =  -T  (5,u  =  -  ~  f (5)  (A. 4. 6) 

5  n  2 

and 


f2q)f'(n) 
5n  f3(C)+f3(n) 


(0  +T?)  + 


f2(pr  (o 

f3(C)+f3(n) 


(e  +1  )  =  0 

n  n 


0(5,5)  =  0^(5, 5)  =  0  (5,5)  =  0 


(A. 4. 7) 


Rather  than  proceed  as  before,  we  shall  now  use  successive  approximations 
as  in  Garabedian’s  book  [10].  We  set 


1 

2 


11  ] 
R 


f  (5) 


f 2(n)f 1  (n) 

f3(5)+f3(n) 


1 

f  (n) 


f2,(i-)f.;(a  l  ,d5dn 

fJ(5)+f  (n) J 


(A.  4. 8) 


0 


n+1 


//(-£ MiAn) 

Rif  (?)+f  (n) 


p  + q  1 

n  f3(e)+f3(n)  nJ 


d5  dr) 


n  >  1 


(A. 4. 9) 


30  30 

where  p  =  -  ,  q  =  - —  and  R  is  the  region  in  the  interior 

n  35  nn  3q 

of  triangle  ABP  in  Figure  IV  .  Then 


' 


);  1 
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t  =  T  +  0 
n  n 

=  T  +  0_  +  (10-9..) 
1  n  1 

and 


at 

n 


+  Pi  +  (Pn-P!) 


(A. 4. 10) 


It  is  well  known  that  for  a  linear  problem  such  as  this  one  that  the 
approximations  do  converge  to  a  solution  to  the  problem.  We  shall  attempt 

to  show  that,  if  e j p-£ I  is  bounded,  that  is,  if  =  0(— )  then 

£ 

jp^-p^l  =  0(e)  n  =  2,3,...  .  (A. 4. 11) 


In  this  case  our  condition 


will,  to  the  order  we  have  been  working  to,  be  given  those  values  of  (5,n) 
such  that 


9T  ^  n 
35  +  pl  =  0 


(A.  4. 12) 


We  obtain,  after  a  short  calculation, 


,c  ^  -  -1  „  )f3(£)+f3(n) 

Pl(?.n)  6f(£)t'n\  2f3(?) 


„,r  ,  1  .  I  f3(p+f3(n) 

1  6f(n)  ^  2f 3 (n) 


f  2 (g ) £  *  (£)  fn  dy _ 

2  ?  f(u)[f3«)+f3(u)] 

(A.  4. 13) 

f2(n)f 1  (ri)  j-n  _ dy 

2  \  £(y)[f3(n)+£3<y)] 


We  wish  to  find  $2*^1  ’  ^2~^1  ’  substitution  of  these  above 

quantities  in  (A. 4. 9).  We  shall  break  this  into  two  parts  P1  ,  Q1 


' ' 


r 


■ 
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II  II 

obtained  from  the  logarithm  terms  in  (A.  4.13j)  and  P  ,  Q  obtained 
from  the  integral  terms.  Again,  after  a  little  calculation  we  find 


P1  = 


1 

6 


f  (y)f'(y) 


f  (?)  [f 3 (5)+f  3  (y )  ] 


£n 


f3(;)+f3(u) 

2f  3 (£ ) 


f  (C)f'  «) 


£(y)[£3(5)+f3(w)] 


Jin  < 


! f3(;)+f3(u) 


2f 3 (u) 


dji 


(A. 4. 14) 


Q1- 


I 

J 


,  ,  1 

f  (n)f’ (n)  on  , 

t(n)+f J(?)  l 

o  o  yn  s 

f(y)[f  (n)+£  (y)] 

2f3(y)  j 

f  (ii)£'(u) 


f  (y) [f3(n)+f3(u)] 


Jin  < 


f3(n)+f 3(y) 


2f 2(u) 


> 


dy 


(A4.15) 


We  claim  that  there  exists  a  positive  constant  K  such  that 


_.I  |  £  K  » 

p  <  — —  ?-n 


■If  £  K  [  ^ 

Q  !  <  -T-  5-n 


To  see  this  we  write,  for  example, 


(A. 4. 16) 


iln  < 


f3(£)+f3(u) 

2f 3 (£) 


=  Jin  < 


1  + 


f3(y)-f3(?) 

2f3U) 


f 3 (y )~f 3  (?) 

2f  3  (? ) 


1+A 


0  <  A  <  1 


(A. 4. 17) 


i  3  i-l  1  1 

Clearly  |f  (?) |  >  — if  we  ensure  that  £  < 


.3/2 


2an 


as  we  do 


here  in  what  follows.  Further 


■ 
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f 3 (u )-f 3 (£)  I 


(l-fealT  cos  Hu)3  -  (1+eaJI  cos  Ilg)3 


(1+eaII  cos  IIy)3^3+(l+ear[  cos  IlC)3^ 


(1+eaII) 3  -  (1-can)3 


2  Crp  3/2 


3  3  3 

6saII  +  2e  a  II 

,-d/2) 


<  K1e 


(A. 4. 18) 


for  >  0  ,  and  e  sufficiently  small.  Noting  also 


that 


f’(x)  = 


-call  sin  IIx 


(1+eaII  cos  Tlx) 


1/2 


(A. 4. 19) 


1  1/2  1/2  1 
and  that  (— )  <  f(x)  <  2  for  all  x  and  all  e  <  -r— -  it  is  not 

l  zaJI 

difficult  to  see  that  the  statement  (A. 4. 16)  is  true  for  some  K  >  0  . 


We  turn  now  to  the  contribution  from  the  second  part  and  find 


II  1 


=  -  f 

(s  J 


n 


f 


f  (z)f'U) 


5  f(y)[£3«)+f3(y)] 


£  (€)£ ' (6) 


6  h  f(y>[f3(y)-£3«)] 


In  < 


Jin 


f3(Q+£3(n) 

f3(C)+f3(y) 


f3(y)+f3(n) 

£3<C)+f3(n) 


dy 


f3q)+f3(y)  \ 

2f3(y) 


(A.  4. 20) 
dy 


while 


.11 


1 

6 


0 


-I  f 

6 


f  (n)f’(n) 


€  f(y)[f3(n)+f3(y)l 


n 


f  (n)f'(n) 


C  f(y)lf3(y)-f3(n)] 


f -3 - 3 


dv 

£J(y)+fJ(y) 


Jin 


2f3(y) 


f3(u)+f3(n)  _  f3(y)+f3(g) 


f3«)+f3(n) 


(A.  4. 21) 
dy 
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The  only  difficulty  is  making  our  estimate  arise  from  the  second  term 
in  each  case.  We  write,  using  (A. 4. 20)  as  an  illustration, 


f3(y)-f3(0 


f3(y)+f3(n)  .  £3(Q+f3(y) 

[f3(S)+f3(n)  2f 3 (y) 


f3(y)-f3(C) 


f3(u)-f 3(g) 

£3(5)+f3(n) 


u 

f3(y)-f3(5) 

2 

• 

-1 

f3(E)+£3(n) 

<1+A>2 j 

f3(y)-f3(C) 

£3(y)-f3(5) 

2 

1 

1 

2f3(y) 

“T 

2f 3  Cu ) 

(m2)2 ; 

(A.  4. 22) 


where  0  <  <  1  .  Then  we  have  this  equal  to 


£3(C)+£3(n) 


2f3(u) 


+  0(e) 


(f3(y)-£3(c)}  +  {f3(y)-f3(n))  +  0(£2) 
2£3(y)  {£3(5)+f3(n)) 


(A.  4.  23) 


It  is  now  a  simple  matter  to  see  that 


* 


[XI , 

Ip  1  i 

2 

£ 

2 

*  i 

K  1 

£-n  |  > 

1  II 1 

Iq  i 

2 

E 

2 

*  i 

K  | 

5-n| 

Taking 

K 

=  max 

:  [ K , K* ]  it 

lp2'pl 

< 

K  £  2 

|S-n| 

k2-q] 

< 

K  £  2 

U-n  |  , 

(A. 4. 24) 


(A.  4. 25) 


The  remainder  of  our  argument  will  now  follow  in  a  standard 
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way.  Let 


Then 


y  =  max 
x,y 


an  sin  nx  f (x) 

f3(x)+f3(y) 


(A. 4. 26) 


1  I 

"Pnl  1  I  /  (|Pn"Pn-ll  +  lqn"Vll)dnl  ’  ?  fixed  ’ 


n+1 


n 

WqJ  ^  ElJ  I  ^  (IWll  +  IVVll)iEl  >  n  flxed  ' 


(A. 4. 27) 


We  then  obtain  for  n  >  1  , 


P  “P 
n 


n-1 1  - 


,0  Nn-2  n 
(2y)  K  e 


;-nl 


:  n— 1 


(n-1)! 


(A. 4. 28) 


and  so 


P  “Pi 
1  rn  ^1 


i  l 

k=2 


CO  | 

<W1>  ■ K  J2(2u)k'2  ek  ^ 


k-1 


(k-1)! 


< 


CO 


(2ye) 


k-1 


(k-1)! 


Ke 

2 


7  [ (2ye) 1 g-n  1  ]n 

4  .  m! 


Ke  r  2ye|?-n 
2  1 


1} 


(A.  4. 29) 


As  we  bound  e  |  £-ri  |  ,  it  follows  that  |p  ^p^| 
D  >  0  ,  and  our  conjecture  that  to  order  e  , 
given  by  the  zeros  of  is  correct. 


<  eD  for  a  finite  constant 
the  zeros  of  t^  are 
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A  simple  calculation  shows  that 


T?  +  Pi  = 

2 

1  an  >■  v  ,  -  \ 

-  2  ”  £  16  (n-5)  Sin  lip  +  0(e)  , 

which  is  in  agreement  with  the  results  of  Chapter  IV. 


■ 


BIBLIOGRAPHY 


Ames,  W.F.,  ’Nonlinear  Partial  Differential  Equations  in  Engineering’, 
Academic  Press,  1965,  New  York,  p.  464. 

Corpson,  E.T. ,  "On  the  Riemann-Gr een  Function’  Archive  for  Rational 
Mechanics  and  Analysis,  Springer  Verlag,  Berlin,  Vol.  I,  1957-58,  p.  324. 

Courant,  R. ,  and  Freidrichs,  K.O. ,  'Supersonic  Flow  and  Wave 
Propagation',  Interscience,  pp.  361,  New  York,  (1948). 

Courant,  R. ,  and  Hilbert,  ’Methods  of  Mathematical  Physics',  Vol.  II, 
Interscience,  1966,  p.  42. 

Cragg,  J.W. ,  'The  Breakdown  of  the  Hodograph  Transformation  for 
Irrotational  Compressible  Fluid  Flow  in  Two  Dimensions,  Proc.  Comb. 

Phil.  Soc. ,  44,  pt.  3  pp.  360-376,  (1947). 

Dienes,  P. ,  'The  Taylor  Series',  Dover,  New  York,  1957. 

Fermi,  E. ,  Pasta,  J.R. ,  Ulam,  S. ,  Los  Alamos,  Report  No.  1940,  May, 

1955. 


Forsyth,  G. ,  and  Wasow,  W. ,  'Finite-Difference  Methods  for  Partial 
Differential  Equations',  Wiley  and  Sons,  New  York,  1960. 

Fox,  P.A. ,  'Perturbation  Theory  of  Wave  Propagation  Based  on  the 
Method  of  Characteristics',  Journal  of  Math,  and  Phys.,  Vol.  34, 
pp.  133-151. 


Garabedian,  P.R. ,  'Partial  Differential  Equations',  Wiley  and  Sons, 
New  York,  1964,  p.  110-117. 


Goursat,  'Differential  Equations',  Vol.  II,  Dover,  New  York,  1959, 
p.  45. 

Lax,  P.D. ,  'Development  of  Solutions  of  Nonlinear  Hyperbolic 
Partial  Differential  Equations',  Journal  of  Mathematical  Physics, 
Vol.  5,  No.  5,  May,  1964,  p.  611. 


Lin,  C.C.,  'On  a  Perturbation  Theory  Based  on  the  Method  of 
Characteristics',  Journal  of  Math,  and  Physics,  Vol.  33,  No.  2, 
p.  117-134,  (1954). 

Miranker,  W.L. ,  'The  Wave  Equation  in  a  Medium  in  Motion' ,  IBM 
Journal  of  Research  and  Development,  Vol.  4,  No.  1,  June  1960. 


Zabusky,  N.J.,  'Exact  Solution  for  the  Vibrations  of  a  Nonlinear 
Continuous  Model  String',  Journal  of  Mathematical  Physics,  Vol.  3, 
No.  5,  Sept.  1962,  p.  1028. 


1 

t  :  •  olVfiT  •  J'J 


103  - 


16.  Zabusky,  N.  J.  ,  and  Kruskal,  M.D.,  ’Interaction  of  Solutions  in  a 
Collisionless  Plasma  and  the  Recurrence  of  Initial  States’,  Physical 
Review  Letters,  Vol.  15,  No.  6,  August,  1965. 

17.  Zabusky,  N. J. ,  and  Kruskal,  M.D. ,  ' Stronboscopic  Perturbation 

Procedure  for  Treating  a  Class  of  Nonlinear  Wave  Equation’ ,  Journal 
of  Mathematical  Physics,  Vol.  5,  No.  2,  Feb.,  1964,  p.  231. 

18.  Zaiser  Janes,  'Nonlinear  Vibrations  of  a  Moving  Threadline’,  . 
Doctoral  Thesis,  U.  of  Delaware,  June,  1964. 


